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HEAT TRANSFER AND LEVITATION OF A SPHERE IN LEIDENFROST BOILING 
by Robert C. Hendricks and Kenneth J. Baumeister 
Lewis Research Center 

SUMMARY 

A hydrodynamic model is postulated for a small, spherical, liquid drop film boiling 
on a cryogenic surface, such as may occur in propellant spillage accidents and in the 
preparation and preservation of some biological species (e.g. , blood and enzymes). The 
momentum and energy equations are solved analytically to predict the velocity and tem- 
perature fields in the vapor film beneath the sphere. This leads to a heat -transfer coef- 
ficient of the form 



+ K 


k 3 p d pg^* 

R 0 M(T d - T s )F(0*) 
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where k is the thermal conductivity of the vapor, p d the density of the sphere, p the 
density of the vapor, X* the modified latent heat, R q the radius of the sphere, p the 
vapor viscosity, T d and T g the temperature of the sphere and the liquid, respectively, 
k a constant, and F(0*) a function which depends on the characteristics of the sphere 
and the supporting fluid. The values of tc and F(0*) and dimensionless forms of the 
heat-transfer solutions are given in the body of the report (see table II). 

This expression for the heat-transfer coefficient combined with a thermal balance 
predicts the time necessary for spheres at room temperature (liquid- liquid film boiling) 
to freeze when placed on a cryogenic liquid (solid- liquid film boiling). The agreement 
between experimental and predicted times is fair. A more extensive freezing analysis 
and thermocoupled spheres appear to be necessary. 

Two solutions for the heat-transfer coefficient between the drop and the supporting 
surface are presented herein, one for a variable vapor-gap thickness and the other for 
constant vapor-gap thickness. These solutions and their associated boundary conditions 
at the liquid- vapor interface revealed two significant facets: (1) Whether the analyst as- 
sumes a constant or variable vapor-gap thickness, the resulting average heat-transfer 
coefficients are not significantly different; and (2) a technique called the integral bound- 
ary constraint can be used to handle undesirable nonconstant boundary conditions. 


INTRODUCTION 


Film-boiling heat transfer between two liquids at different temperatures has been 
under extensive investigation recently in consideration of propellant spillage accidents 
(ref. 1). If liquid propellants spill accidently during test-stand or launching operations, 
a catastrophic explosion resulting from the detonation of the fuel and oxidizer is possible. 
In particular, a number of experiments (ref. 1) have been performed involving the mixing 
of small kerosene spheres in liquid oxygen. When small kerosene spheres are dispersed 
in a matrix of liquid oxygen, the sensible heat of the kerosene spheres vaporizes the 
liquid oxygen, thereby forming a film of oxygen vapor around the sphere. This vapor film 
has a significant damping effect on the mixture's explosive potential, since gaseous oxy- 
gen contacting the liquid fuel has a much lower explosive potential than concentrated 
liquid oxygen. 

The liquid-liquid film-boiling phenomenon is also of interest purely from stability 
considerations. In this case, a liquid sphere with a higher specific gravity floats upon a 
a cryogenic liquid of lower specific gravity while film boiling is taking place. The sphere 
usually freezes and continues to float (solid- liquid film boiling); suddenly at the transition 
from film boiling to nucleate boiling, the vapor blanket surrounding the sphere disappears. 
Then the sphere sinks beneath the liquid and falls to the bottom of the pool. 

This report considers two aspects of this liquid-liquid film-boiling problem. First, 
a hydrodynamic model is postulated for a small, spherical, liquid drop film boiling on a 
cryogenic surface. The conservation equations of momentum and energy are solved anal- 
ytically to predict the velocity and temperature fields, and a simple theoretical correla- 
tion for the heat-transfer coefficient is presented. 

Actually, two solutions are presented, one for the case of variable vapor-gap thick- 
ness, which is limited to 0* < tt/2 (see fig. 1); and the other more general solution for 
an assumed constant vapor- gap thickness, which is valid for 0* < 77. The angle 0* is 
the interface separation angle, that is, the angle from the stagnation point to where the 
sphere and the interface separate. This angle is discussed in detail in the next section. 
Two completely different similarity transforms are used to obtain the different solutions; 
and, to attain an exact solution to the variable -vapor -gap problem, a technique of condi- 
tioning the boundaries called the integral boundary constraint is used. 

Both solutions, however, are shown to yield expressions for the heat-transfer coef- 
ficient which differ only by a constant which is much smaller than can be detected by our 
experiment. The agreement between the two solutions gives some credence to the almost 
universal assumption of constant vapor-gap thickness in solving film-boiling problems 
(refs. 2 to 4). 

Secondly, a heat balance is performed on the sphere, which leads to a correlation 
equation for predicting the time required for a liquid sphere to cool from room tempera- 
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ture to its freezing point and then to the temperature at which boiling ceases; liquid sub- 
cooling was neglected. Theory is checked with experiment by placing water spheres on 
liquid nitrogen and comparing the experimental and theoretically predicted times for the 
liquid sphere to freeze. Fair agreement was found. 

A motion-picture supplement C-267 has been prepared and is available on loan. The 
film supplement demonstrates and discusses the phenomena associated with levitated 
spheres in Leidenfrost boiling. A request card and a description of the film are included 
at the back of this report. 


BASIC MODEL AND EQUATIONS 

Consider the film-boiling model depicted in figure 1. This model appears to fit the 



physical case of a solid or liquid sphere "floating” on a second liquid (solid-liquid and 
liquid-liquid film boiling); for example, drops of water floating on a sea of liquid nitro- 
gen. There is, of course, no reason to limit the model to liquid spheres or cryogenic 
fluids. Slush and small solid materials could also be utilized on any liquid provided the 
proper temperature difference and buoyancy criteria are maintained. 

The model used in this study applies to two cases with spherical geometries: 

(1) In the first case , the vapor-gap thickness for a floating sphere is permitted to 
vary, although the subsequent solution is limited to 8* < n/2. Such a solution seems to 
model the physical process quite well. However, some difficulties are encountered, and 
the limitation to 9* < u/2 is too constrained. It will be shown that the average heat- 
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transfer coefficient for the variable-vapor-gap case vindicates the assumption of constant 
vapor-gap thickness; however, the authors sought another solution to the problem, which 
is applicable to 9* < v . 

(2) The second case , which assumes a constant vapor-gap thickness, sacrifices some 
of'the physics (e.g. , the variable gap thickness) in return for a solution which is appli- 
cable for 9* < 77. 

The consideration of these two cases requires different similarity transforms, as dis- 
cussed in the section Similarity Transforms. The model may be applied either to the 
case where the sensible heat of the floating sphere evaporates the encompassing liquid or 
to the case where the sensible heat of the encompassing liquid evaporates the floating 
sphere (e.g. , liquid-nitrogen sphere on water). For simplicity, it is assumed throughout 
the remainder of this report that the supporting fluid evaporates at the interface. Phys- 
ically, for the experimental results considered herein, the heat passes from the warmer 
spherical drop and evaporates the fluid material beneath it, thereby forming a supporting 
vapor layer under the sphere. The vapor layer is, in turn, supported by the interface. 
Surface tension, density ratio, buoyancy, surface curvature, and the "nonwetting" char- 
acter of the interface play an important role in determining whether the interface is 
strong enough to support the sphere. During the evaporative, or film-boiling, process 
the temperature of the sphere decreases until it coincides with the Leidenfrost tempera- 
ture of the liquid. At this time, film boiling ceases, and transition boiling quickly fol- 
lowed by nucleate boiling begins. The vapor layer beneath the sphere no longer exists; 
rather, small nucleate bubbles form in the sites on the solid (liquid) sphere. Thus, be- 
cause there is no longer a supporting vapor gap, the sphere sinks beneath the surface and 
falls to the bottom. 

The following assumptions are made in developing the model: 

(1) Evaporation and internal circulation of the liquid sphere are considered to be 
small. 

(2) The model has complete symmetry with respect to the angular coordinate <f>. 

(3) The flow of vapor is considered to be laminar and incompressible; and the inertia 
terms in the Navier-Stokes equations are neglected. Justification for this assumption can 
be found in reference 2. 

(4) Radiation is assumed to be negligible. 

(5) The velocity and temperature profiles are assumed to be in steady state. 

(6) At any instant of time, the sphere is at an average temperature T^, and the 
evaporating liquid is assumed to be at the saturation temperature T_. The properties of 

fa 

the flow field are evaluated at the film temperature 


and are considered to be constant. This assumption has worked quite well (e.g. , ref. 5). 
(All symbols are defined in appendix A.) 

(7) The convective terms in the energy equation, as treated in appendix F, effect a 
correction factor on the conduction term. The correction factor is taken to be 1 (i.e. , 
conduction dominates) to make the problem tractable. Such an assumption is apparently' 
a good one based on the work in references 2, 5, and 6. These authors indicate the major 
mode of heat transport to be conduction, as further verified in appendix F. Thus, it is 
assumed that 


^ 9T <y 3T «__k_ ± L2 3T\ 
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(8) Heat transport within the sphere is by conduction alone. With cooling at the lower 
surface, there should be no instabilities or ensuing cellular motion. 

(9) The supporting liquid is at the saturation temperature. Thus, with assumption 7, 
all the heat reaching the liquid produces vapor. 

(10) Subcooling of the liquid sphere was considered but omitted because of the presence 
of the freezing crystalline ice shell. 


Governing Equations 

The governing equations are as follows (ref. 7): 
Momentum: 


o =lf£ + dv 2 v - . i. Hi - 2 . v„ cot «1 

P 3 r V r 2 2 30 2 


-g c 1 9P / 2 2 ^r 

0= — E i — + v V z V a +- 


p r 30 


0 


2 30 2 . 2 a . 

r r sin 0> 


(3) 

(4) 


where 
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1 


(5) 



_3_ 

3r 



r 2 sin 



Energy: 



Continuity: 


(6) 


0 = -L JL (r 2 V ) + - (Vn sin 9) 

2 3r r sin 9 d 9 

r 


(7) 


The boundary conditions on the velocity and temperature fields are given in table I, 
where 9* is defined in figure 1, and 5 is the vapor- gap thickness . 


TABLE I. - BOUNDARY CONDITIONS FOR FLOATING SPHERE 


Area 

Surface 

Boundary condition 

Equation number 

Sphere surface 

r = R o 

eie* 

v r = 0 

v 0 = o 

T = T d - e(T d - T g ) cos 0 

(8) 

Liquid- vapor 
interface 

r = R 0 + 5 

e g e* 

V r = V r (6) 

v e = o 
T = T S 

(9) 

Stagnation 

region 

R 0 £r - R o +6 
0 = 0 

II 

o 

(10) 

Separation 

region 

r = R 0 
0 =* 0* 

p.p 0 

(11) 


The tangential Vq and radial V r vapor velocities are both zero at the surface of 
the sphere (eq. (8)), while only the tangential velocity Vq is zero at the surface of the 
cryogenic supporting liquid (eq. (9)). Because of evaporation, a radial velocity V r (5) 
exists at the surface of the cryogenic fluid. 
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The cryogenic fluid is assumed to be at the saturation temperature T , as shown in 

s 

equation (9) . The boundary condition at the surface of the sphere is more complicated 
and, in fact, an unknown. In this report, the general case of a surface temperature of 
the form 


T = T d - e (T d - T s ) cos 9 (11a) 

is investigated. Specifically, two cases are considered: 

(i) The sphere is assumed to be at an average temperature, a condition discussed 
more fully in the section Integral average constraint (p. 9); that is, using equation (11) 
and integrating over the surface of the sphere, 

T d =i/ TdA = T d (llb) 

(ii) The temperature at the surface of the sphere follows a cosine variation with some 
assumed value of e (a positive constant which is less than 1). This is investigated for 
the first case only (i.e. , when 6* < u/2). 


Constraints 


The boundary conditions are still incomplete at this point since V r (6), 6, and 9* are 
unknowns. Hence, three additional mathematical constraints are necessary to make the 
problem tractable . A fourth mathematical constraint is required in the case of the vari- 
able vapor-gap thickness to force a tractable solution. 

Static force balance (neglecting forces internal to sphere) . - One additional con- 
straint requires the sphere to be in static equilibrium. Small vibrations of the sphere 
are neglected. The static equilibrium condition requires that the weight of the sphere be 
balanced by the shear and pressure forces acting beneath the sphere. These forces are 
depicted in figure 2. Summing the forces acting in the vertical direction gives 


°= W d- 


/*9* 

*/n t/n 


(P cos 9 + t Q sin 9 - r rr cos 9) 


r=R, 


R* sin 9 d 9 d$ 



P Q cos 9 R^ sin 9 d9 d<b + 7rR^P Q 


(12) 
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The last two terms in equation (12) represent ambient pressure forces acting on the 
sphere. 



Figure 2. - Schematic model of pressure and shear forces acting 
on sphere. 


Interface energy balance . - The second constraint necessary for the solution of equa- 
tions (3) to (7) is the interface energy balance. Because the supporting fluid has been 
assumed to be at the saturation temperature (assumption 9), all the heat leaving the 
sphere produces evaporation of the fluid material. Mathematically, this constraint is 
expressed as 


-P AV r (R o +6,9) = -k — 

ox 


r= R 0 *> 


(13) 


where A is the latent heat of vaporization and 



r= R o+ 6 


is the conduction heat flux 


to the boundary of the supporting liquid. 

Free-surface pressure head . - The third constraint is that the forces supporting the 
sphere must be transmitted and balanced by the supporting fluid interface. Here the free- 
surface heat Z as shown in figure 2 (see also fig. 8), must be determined from a bal- 
ance of forces acting on the liquid- vapor surface membrane. The determination of Z Q 
such that all forces are balanced and the sphere is supported also uniquely prescribes a 
value to 9 *. The solution for Z Q is discussed in appendix H. 
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Integral average constraint . - A fourth constraint which the authors found to be a 
very valuable tool in solving problems with nonconstant boundary conditions is constrain- 
ing the problem to satisfy a suitably averaged condition at the boundary. While the tech- 
nique has general applicability, the constraint was used for the following two cases. 
First, for the variable-vapor-gap case, it was found that in order to make the problem 
rigorous, the radial velocity at the liquid-vapor interface had to be a suitable average 
value if the surface temperature was assumed to be constant. Employing the integral 
average constraint yields 


V fi =- / V (5 , 9) dA (13a) 

0 A *A r 

In similar manner, the average surface temperature is defined by using equation (11) and 
the integral average constraint 


T 



T (1,0) dAj 


(13b) 


where A^ is not necessarily equal to A, a choice which is available to the analyst. For 
example, A^ could represent only the area influenced by the vapor (0^0 s 0*), in which 
case 



/T dA 
A 


2tt 


/<v 


€ AT cos 0)R^ sin 9 d 9 


2 ? t 



sin 9 d 9 


T , = T , - 1 — (1 + cos 0*) 
d d o 


(13c) 


If, however, the analyst chooses A = Aj (0 ^ 0 — v ) , then = T^. 

The solution of the preceding set of differential equations is given in the next section. 
Figure 3 presents a solution flow chart of the problem, which may be used as a study aid. 
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Flow 

equations 


Similarity 

transforms 


Energy equations 


Nondimensionalization 

: 


n 



Radial velocity 
v A « C 4 A 3 cos 0 


| Interface velocity 
distribution 


Ordinary 

flow 

equations 


"C 4 " 

value 



Static 


force 


balance 


Nondimensionalization 


V A H 


Interface 

energy -{A}^{h] 
balance 


Determination of 
free- surface 
height Z 0 


If sphere will not float, 
see ref. 10 for submerged- 
sphere analysis 


Figure 3. - Solution flow chart 


SOLUTION 

Momentum Equations 

The governing equations are made nondimensional to generalize the solution by sel- 
ecting the following parameters: 


C = — 

R o 

V r 

v r = — 

s u* 


Vn = 


u* = 


R 


(14) 


(15) 


(16) 


(17) 
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( 18 ) 


pu*R 


Re = 2 h l 

M 


P =■ 


p u* 


(19) 


Substituting these parameters into the momentum equations gives 


0 = -— + V 2 V;»- — V. - — - — v n cot 0 

5 ,2 * 


a? 


^2 90 ^2 


( 20 ) 


0 = -_ 

c 90 


1 0J) + v 2 v. + ±- 


2 9v C 


v„ 


<? 00 C 2 Sin 2 0 


( 21 ) 


where 


V 2 = J_ JL 0 \ i 


C 


2 9 c; 


sin 8 


L_ JL (sin 0 ±\ 
... „ ae V del 


( 22 ) 


Introducting the stream function for spherical coordinates (ref. 7, p. 131) gives for 
the velocity distribution 


V r = - 


dip 


r 2 slnS 36 


(23) 


y - 1 

0 : ~ T 

r sin 0 9r 


(24) 


or, in terms of the dimensionless velocities, 


v c = - 


i MM) 


C 2 sin 0 cci 


(25) 


Vn =• 


J W) 

£ sin 0 9$ 


(26) 
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Substituting equations (25) and (26) into the momentum equations (eqs. (21) and (22)) 
and combining to eliminate the pressure terms give 

E 4 (!//) = 0 (27) 


where 


E 


2 


8 2 [ sin e J_ / 1 

^2 ^2 d8 \sin 0 



(28) 


Similarity Transforms 

The governing stream function differential equation (eq. (27)) is now reduced to an 
ordinary differential equation by means of similarity transforms. For the variable- 
vapor-gap-thickness case, the transform is 

Tj: i/yC,0) - Fj(0 sin 2 6 (29) 

which is restricted to those values of 6 less than v/2. For the more general but less 
physically realistic case of constant vapor-gap thickness, the transform is 

T 2 : tf/ 2 (C,0) = F 2 (0(1 - cos 6) (30) 

For small spheres where the density ratio between the sphere and the supporting 
fluid is near 1, the value of 0* is less than u/2. In this case, either transform gives 
approximately the same value of the heat-transfer coefficient. At present, which solution 
more nearly approximates the actual flow has not been determined. 

For heavy or large-diameter spheres, however, 6* is greater than tt/2. In this 
case, only the second transformation can be used. 


Heat -Transfer Coefficient 

The details of the solution of the momentum equation and the energy equation are 
long and involved and are performed in the appendixes. 

By using the similarity transform T^ (eq. (29)), a solution of the governing equations 
can be attained. Appendix B gives the solution for the heat-transfer coefficient when 
0* < 77/2 for two cases: (1) constant surface temperature and (b) strong surface- 
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temperature variation with angular position 0. 

For d* beyond n/2 (and less than 77 ) , the similarity transform T 2 (eq. (30)) ap- . 
plies and a solution for the heat-transfer coefficient is presented in appendix C. 

The final expressions for the heat -transfer coefficients are presented in table II. 

In the expression for the heat-transfer coefficient, the effect of convection on the 
energy equation was neglected. Appendix G shows that the effect of convection will be of 
the order of 10 percent. At the expense of accuracy and for increased simplification, 
however, the closed-form conduction solution was used in the analysis, and A was re- 
placed by A* (eq. (G16)). 


TABLE II. - SUMMARY OF SOLUTIONS TO FLOATING SPHERE PROBLEM 


Transform 

Boundary condition 
on surface 
temperature 

Nusselt number 

Constant, 

K 

F(0*) 

Equation 

number 

T 2 : F 2 (0(1 - COS 9), 
9* < n 

Constant at T d 

Nu d = 1 + 
f(9*) = 

2 Ra'/ Pd \ 1 

9 \p t - pj f (9*) 

1-2 cos 0* + cos^ 0 

1/4 

* 

(!f 

f(9*) 

(C25a) 

Ty F j(5) sin 2 9) , 

e* < 77/2 

Constant at 

Nu d = 1 + 

G(9*) = - 
1 

1 

Ra' / Pd \ 

ec(e*)\p l - P j 

g(0*) 

+ cos 9 * 

- 1 cos 0* + i cc 
2 2 

1 + cos 9* 

1/4 
)s 3 9* 

(if 

G(9*) 

(B45a) 

Strong variation 
with angular 
position 0 

Nu d = l + 

where cc 
trary c 

Ra' / p d \ cos 0 ref 

1/4 

/2 

/ cos 9 ref\ V4 

g(0*) 

(B66a) 

3g(**)\ft - pj <4 

s is some arbi 

onstant, usually near 1 

l 3e ! ) 


COMPARISON OF EXPERIMENT AND THEORY 

A simple test of the analytical model is performed using water and liquid nitrogen. 
Liquid nitrogen represents the supporting surface, called the evaporative fluid material, 
shown in figure 1. When small quantities of water are placed on the surface, the water 
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forms a spherical drop because of the action of surface tension. At the expense of the 
internal energy of the water spheroid, liquid nitrogen is vaporized beneath the sphere. 

This vapor forms the supporting cushion for the sphere, which appears to "float" on the 
much colder nitrogen surface (liquid- liquid film boiling). 

During the experiment, a slight amount of dye coloring was added to the water. Dur- 
ing the freezing process, the dye color changes and continues to change until the entire 
drop is frozen. The frozen drop is now a hard sphere "floating" on a sea of liquid nitro- 
gen (solid-liquid film boiling). Its free-floating state is sharply terminated at the onset 
of nucleate boiling, and the sphere falls beneath the surface. 

The size of the water sphere used can be determined by measuring the residue (the 
hard sphere of ice) or by using a pipette. Both techniques were used. A stop clock was 
used to determine the elapsed time between placing the water on the nitrogen and the color 
change in the sphere, time for the sphere to be completely frozen, and the time for the 
sphere to fall to the bottom. 

In performing the experiment, the water sphere fell about half a centimeter onto the 
liquid nitrogen. As film boiling must be initiated (negligible radiation) through contact of 
the participating elements, in this case the two fluids, a thin layer of ice could form be- 
neath the sphere. If this hypothesis is correct, the ice layer would grow very slowly until 
the internal temperature of the liquid were sufficiently low, at which time the ice layer 
would advance quite rapidly. It seems there would be little convection, in light of the 
stability provided by the temperature gradient. There is, however, a natural inverse 
gradient at 4° C and liquid subcooling is certainly a factor. The inverse gradient could 
enhance the freezing rate because of the density inversion occurring at this temperature; 
conversely, the subcooling would degrade the freezing rate. These effects could be self- 
cancelling. 

To use the experimentally measured sphere freezing times as a check on the theoret- 
ically determined heat-transfer coefficients, it is necessary to determine the temperature 
of the sphere as a function of time during freezing. Two models are used. The first, a 
simplified model, assumes that no thermal gradients exist in the sphere. This is analo- 
gous to the classical Newtonian cooling problem. The second model, called the pseudo- 
steady-state model, takes into account thermal gradients in the freezing sphere (no sub- 
cooling). The calculated time-temperature histories in the freezing sphere for both 
models are shown in figure 4. The derivation for the pseudo-steady-state model (lower 
pair of curves in fig. 4) is given in appendix I. Only the simplified model for sphere tem- 
perature is presented in the main text. 
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Figure 4. - Calculated surface temperature for a 61-centimeter-diameter water sphere on liquid ni- 
trogen as function of time. 


In the simplified model calculations, the surface of the water sphere is assumed to 
be at 0° C and to remain at that temperature until the remainder of the sphere reaches 
0° C. Therefore, neglecting the mass of the small, initial, free ice layer, the sensible 
cooling time can be found from an energy balance on the sphere 


< hA >d< T f - T s>‘sen - m ( C p)d (T d,o - T,) 


(31) 


Solving for the time to remove the sensible heat t gen gives 

”(S) d < T d.o - T f> 

sen ' (hA) d (T f - T g ) 


(32) 


For a 0. 61- centimeter-diameter water sphere on nitrogen, the value of t is 4.2 sec- 

S6u 

onds, as shown in figure 4 (see eq. (J21), appendix J, for the numerical calculation). 

The time required to freeze the entire sphere isothermally can be calculated in the 
same manner, only C AT is replaced by the latent heat of fusion y, that is, 

Ir 


m ice y 


ice 


(hA) d (T f 


T s> 


(33) 


For the 0. 61- centimeter- diameter water sphere on nitrogen, the value of t^ ce is 12. 5 
seconds. This is represented by the horizontal section of the upper curves in figure 4. 
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Once the entire sphere has turned to ice, the solid sphere loses its sensible energy 
and cools to the temperature of the supporting liquid. The sensible cooling time can be 
found from an energy balance on the sphere, of the form 


/ \ S ( T d " T s) 

h d A (T d -T s )=-( P C p ) (ce V. 


3t 


(34) 


where 


r = t - (t 


sen 


+ t ice ) 


(35) 


and at r = 0, the surface temperature T(0) is taken as the freezing temperature T^. 
The solution of equation (34) becomes (for an average value of ^/(^p). ^ ’ 




(36) 


or solving for the time r, 


T = 




(37) 


Thus, the time for a sphere to freeze to some temperature beneath its freezing point is 


t 


= t 


sen 


+ 1 . 


ice 




(38) 


From this equation, the temperature falloff shown by the upper curves in figure 4 is pre- 
dicted. 

As mentioned earlier, as the free-floating sphere cools and its outside surface tem- 
perature approaches the saturation temperature, film boiling ceases and nucleate boiling 
begins. At the onset of nucleate boiling, the sphere is wetted and falls beneath the sur- 
face. The transition temperature at which nucleate boiling occurs is called the Leidenfrost 
temperature. Assuming, for a moment, that the Leidenfrost temperature is accurately 
known, the important temperature in predicting the transition temperature is then the 
outside surface temperature of the sphere T Q . It can be shown that 
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( 39 ) 


T < T h 
o d 


That is, at a given time, the outside temperature of the sphere T Q will always be less 
than the temperature T d calculated from the simplified analysis. This is because a 
temperature gradient always exists throughout the sphere. Consequently, the more exact 
pseudo -steady- state analysis given in appendixes I and J predicts the lower curves in 
figure 4. In either case, the sphere is assumed to be smooth. 

The Leidenfrost temperature for water on nitrogen was estimated to be 126.0 K (see 
appendix K for estimation). This is illustrated by the line labeled T^ eg j. in figure 4. 

In determining T^ consideration was made for surface material and roughness (ice 
protrusions). The surface temperature of the 0.61-centimeter-diameter water sphere 
falls to this value in 29.2 seconds, at which time transition boiling begins and the sphere 
sinks. The theory is in reasonable agreement with the data of table III (see also fig. 4). 

Additional calculations for other key times associated with the freezing of the sphere 
are given in appendix J. 


TABLE III. - FREEZING DATA ON SMALL SPHERICAL WATER DROPS FLOATING ON 
LIQUID NITROGEN AT AMBIENT CONDITIONS 
[Initial sphere temperature, 298 K.] 


Volume of sphere, 

V 

cm 3 

Radius of sphere, 

R o’ 

cm 

Time for dye 
change to 
become 
visible, 
sec 

Time for sphere 
to be completely 
frozen, 
sec 

Time for sphere 
to fall to 
bottom , 
sec 

a 

0.09 

0.278 

7 

16 

27 


. 113 

a . 3 

8 

17 

23 

a 

.075 

.26 

9 

17 

28 

a 

. 11 

.295 

9 

18 

27 


Measured value. 


CONCLUSIONS 

Theoretical expressions for the heat-transfer coefficient between a liquid sphere 
floating in film boiling and the supporting cryogenic fluid were derived from the funda- 
mental equations of momentum and energy. The expression for the heat-transfer coeffi- 
cient and the pseudo-steady-state freezing analysis lead to an accurate prediction of the 
time necessary for a sphere at room temperature to freeze when placed on the cryogenic 
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fluid. Further analysis of the freezing process (which is quite complex) and the use of 
instrumented spheres are required. 

The heat-transfer coefficient to a sphere floating on the surface in film boiling is 
shown to be equal to 



k 3 gp d pA* 
p(T d - T s )F(0*) 


1/4 


where k is thermal conductivity, R Q radius of sphere, k a constant (see table II), g 
acceleration of local gravity, p d density of sphere, p density of vapor, A* modified 
latent heat of vaporization, p viscosity, T d sphere temperature, T g saturation temper- 
ature of liquid, F(0*) represents f (0*) or G(0*) for 6* < ir and 0*<7t/2, respec- 
tively. These functions and nondimensional forms of the heat-transfer solution are given 
in table II. As seen in this equation, there are two distinct contributing factors to the 
heat-transfer coefficient. The first term represents the heat conduction from a sphere 
to a stagnant fluid, while the second term represents a contribution similar to that for 
parallel- surface Leidenfrost boiling. 

The stagnant conduction will dominate whenever the radius of the sphere R Q is such 

that 


R 0 < K 


4/3 k ^ T d T s> 

_p d pgAF(0*)J 


1/3 


Two separate and distinct transitions are seen in the film boiling of a room- 
temperature liquid sphere on a cryogenic surface. First, the sphere cools from room 
temperature and freezes. For water drops containing a blue dye, the freezing process 
becomes visible to the eye by a very noticeable change in the dye color. Secondly, when 
the sphere temperature cools into the nucleate boiling region, the supporting vapor layer 
is destroyed and the sphere falls beneath the surface. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 7, 1970, 

129-01. 
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APPENDIX A 


SYMBOLS 


A 

surface area of sphere 

Sc 

gravitational constant 

‘c.’Vc 

Ci r) O 

constants 

constants 


in Newton's law of 
motion 

ill 


g(0*) 

function defined by 

C P 

specific heat 

eq. (D26) 

( C p) d 

specific heat of sphere 

h d 

heat-transfer coefficient 

C R 

constant defined by 


to sphere 

eq. (120) 

h i 

heat-transfer coefficient 

D 

differential operator 


at interface 

E 2 

operator defined by 
eq. (28) 

I l,2, 3,4 

nondimens ional compo- 
nents of force 

F 

function defined by 

2, 3, 4 

components of force 


eq. (B2) 

K 

index 

f 2 <0 

function defined by 

k 

thermal conductivity 


eq. (C2) 

L 2 

characteristic length 

F(0*) 

used to represent either 
f(0*) or G(0*) for 

m 

mass of water sphere 


9* < 77 or 9* < 7t/2 , 
respectively 

m c 

arbitrary constants, 
m = 1,2,3,... 

F a 

surface-tension force, 

m ice 

mass of ice formed 


eq. (HI) 

Nu d 

Nusselt number at 

F b 

buoyancy force, eq. (H2) 


sphere-vapor inter- 
face, h.R /k 

f(0*) 

G(0*) 

function defined by 
eq. (C30) 

g(0*)/(l + cos 9*) 

Nu q 

Nu. 

a o 

Nusselt number at sur- 
face of sphere 

Nusselt number at vapor 

Gr 

Grashoff number, 
eq. (B69) as 

p(p t - p)Hpg/ M 2 

1 

liquid interface, 

h iV k 


N * 

membrane force acting 

g 

acceleration of local 
gravity 

in meridian plane, 
fig. 9 
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N e 

membrane force acting in par- 

^sen 

time necessary for sensible 

allel circle plane, fig. 9 

heat to be removed from 

p 

pressure 


sphere 

P o 

atmospheric pressure 

u* 

dimensionless velocity, i>/R q 

Pr 

Prandtl number 

V d 

volume of sphere 

P 

dimensionless pressure, 

V r 

radial velocity component 


P/(pu* 2 /g c ) 

V 5 

radial velocity at liquid- vapor 

Pa 

dimensionless atmospheric 


interface 


pressure, P 0 /(pu* 2 /g c ) 

V 6 

tangential velocity component 

R ^ y R^ 

Gaussian radii of curvature 

v ? 

dimensionless velocity in radial 

R 3 

radius defined in fig. 9 

direction, V r /u* 

Ra' 

modified Rayleigh number, 

V A 

velocity at interface 


eq. (B70) 

v e 

dimensionless velocity, in tan- 

Re 

Reynolds number, pu*R Q /p 


gential direction, V^/u* 

R o 

radius of sphere 

w d 

weight of sphere 

r 

radial coordinate 

w* 

nondimens ionalizing factor, 

o 




pu*R^/g 

S 

property group defined by 


U L. 


eq. (B30) 

y f 

surface force, fig. 9 

T 

temperature 

y 

dummy variable 

T d 

sphere temperature 

z 

coordinate, fig. 9 

Nref 

reference temperature, 
eq. (F10) 

Z o 

z 

initial pressure head, fig. 8 
dimensionless length, Z/L 

T f 

film temperature, (T d + T g )/2 

0 

dimensionless thickness, 

T o 

surface temperature during 


1 - (r f /R 0 ) 


freezing, see appendix I 

r 

dimensionless temperature 

T s 

saturation temperature of liquid 


defined by eq. (15) 

t 

time 

v 

latent heat of fusion 

t. 

time for heat of fusion to be 

A 

dimensionless gap thickness, 

ice 

removed from sphere 


6/R 0 

^obs 

observed time for sphere to 

a a 

approximation for A, trans- 

freeze 


form Tj 
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6 

e 


e 

0 

0 c 

e 

e t 

ref 


9 * 

K 

X 

X* 


A 


M 

v 

I 

^steel 

^ice 

P 


P d 


similarity form for A, trans- 
form Tj 

vapor -gap thickness 

constant, 0 < e < 1, fraction 
of temperature difference 
(T d - T s ) 

dimensionless radial coordinate, 
r / R 0 

dimensionless temperature, 

(T - T„)/(T S - T d ) 


dimensionless temperature, 

<T -V/N ref -T s ] 

angular coordinate 

angular position for evaluating 

( T a) ref - ««• < F8 > 

defined in fig. 1 
constant, see table II 
latent heat of vaporization 
modified latent heat of vapori- 


c p< T d 



zation, a{l+— 
2 


convection correction term, 
eq. (Gil) 

viscosity 

kinematic viscosity 

I evaluated for stainless -steel 
properties 

£ evaluated for ice properties 
density of vapor 
density of sphere 
density of liquid 


a 

T 

T r 9 ’ T rr 
<$ 

<p 

* 


density of material 

surface tension 

dummy variable 

shear stresses, see fig. 2 

dimensionless shear stresses 

angular coordinate defined as 
used, see figs. 2 and 9 

dimensionless function, 

1 + A = 1 + (6/R 0 ) 

stream function 

w d /[87rw*sg(e*)] 


w angular coordinate, see fig. 8 

Subscripts: 

A approximation for A 

amb ambient 

d sphere 

f freezing 

fb film boiling 

ice frozen material 

L Leidenfrost 

l liquid 

o initial 

obs observed 

S similarity form for A 

s saturation 

r radial direction 

ref reference 

tot total 

v vapor 
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Superscripts: 

— average value 

t derivative with respect to 

independent variable 
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APPENDIX B 


SOLUTION OF HEAT-TRANSFER COEFFICIENT FOR FIRST 
SIMILARITY TRANSFORM Tj (WHERE 0* < tt/2) 

The momentum and energy equations were presented in the body of this report along 
with the final solutions. The details of the analyses are presented in this appendix. 


MOMENTUM EQUATIONS 

The governing momentum equations for the flow in the r and 0 directions beneath 
the drop (see fig. 1) are given by equations (3) and (4) in the body of this report. The 
momentum equations were made nondimensional (eqs. (20) and (21)) and then combined 
into a single fourth-order equation by the use of the stream function 

E V = 0 (27) 


where 


g2 = 3^ + sin 6 8 / 1 

-^2 ^ 2 30 \sin 0 

and the stream function \ fy is defined in equations (23) and (24). 

In order to solve equation (27), it must first be converted into a ordinary differential 
equation. The first of two possible similarity transforms 

^(5,0) = F 1 (C)sin 2 0 (29) 

is used to convert equation (27) into an ordinary differential equation. Another possible 
similarity transform is discussed in appendix C. The transform given by equation (29) 
has also been used for solving the problem of flow around a sphere (ref. 7, p. 132). 
Substituting the similarity transform given by equation (29) into the governing equation (27) 
yields 

,2 r^)( D rj) F > (?, = 0 <B1) 



-) 

30 / 


(28) 
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Equation (Bl) is a linear fourth-order homogeneous equation, the solution of which is 


Fi«) =Y + C 2 ? + C 3 ? 2 + C 4 C 4 


(B2) 


The expressions for the dimensionless velocities and Vq are given in terms of 
the stream function i j/ by equations (25) and (26) in the body of this report. Substituting 
equation (B2) into equations (25) and (26) gives 


v r . 


-2Fj cos 6 


(B3) 


v 0 = 


F'j sin 9 


(B4) 


The constants in the previous equations must be determined from the boundary con- 
ditions (eqs. (8) and (9), table I in the body of the report; note that eq. (10) is satisfied 
by (B4)). In terms of the dimensionless variables, these boundary conditions become 


where 


v 0 (l,0) = O 

(B5) 

V 0 (l + A, 9) = 0 

(B6) 

v ? (M) = o 

(B7) 

Vj»( 1 + A , 9 ) = V A 

(B8) 

A,i 

(B9) 


Here, A represents the dimensionless vapor-gap thickness concentric with the sphere. 

Applying these conditions allows the four constants to be expressed in terms of one 
constant of integration, as follows: 


2C, 


c 1 = - 


1 + cp + 4cp* 


[ 


(p 2 + 2 cp 3 + 2 cp 4 


1 


(BIO) 
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where 


Cg = [~1 + cp + cp^ + 6 (<£> 3 + <p 4 ) 

1 + cp + 4 cp^ 


(BH) 


C 3 = - 


1 + cp + 4 cp* 


3+3 cp + 8 (cp^ + <p 3 + cp 4 ) 


(B12) 


cp = 1 + A 


(B13) 


The dimensionless gap thickness A is assumed at this time to be an unknown con- 
stant. However, as is shown later (e.g. , see eq. (B27)), this requires that the interface 
energy balance be satisfied on an average over the heating surface, rather than at every 
point. 

The constant C 4 can be found by applying the static force balance constraint 
(eq. (12)). First, however, the pressure distribution in the vapor gap must be expressed 
in terms of the unknown constants. Substituting equations (B3) and (B4) into equations (20) 
and (21) and solving for the pressure distribution give 


/ 4F , 2F’\ „ 

p(C,0) = - cos ejFY' + — i] + 6 C 

^3 j»2 


(B14) 


which upon substituting in the value of Fj gives 


2c 

p(C,0) = - cos 6 (20C 4 ? + — 1) + 3 C 


(B15) 


where C 2 is related to C 4 by equation (B 11) . 

Now, the constants C 4 and 3 C are determined from the static force balance in 
appendix D. This gives 


C 4 = 


'JV\ 


1 + <p + 4 cp* 


[l + cp +cp l 2 + 6((p 3 + <p 4 )J g(0*) 


(B16) 


l where 
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l 


(D26) 


g(0*) = 1 - 3 cos 0* +- cos 3 0* 
2 2 


and 


3 C = p a + cos 0* (20C 4 + 2C 2 ) (B17) 

2 2 

where w* is defined by equation (B25) as w* = pu* R Q /g c . 

Therefore, utilizing C 4 in Cj, C 2 , and C 3 , the velocity and pressure distributions 
are known relations in cp where cp = 1 + A. 

The energy equation is considered in the next section. Afterwards, the solutions to 
the energy and momentum equations are combined in the interface energy balance to ob- 
tain a solution for the heat-transfer coefficient. 

ENERGY EQUATION 

The energy equation is solved in appendix F for both constant surface temperature 
T = T^ at r = R Q (eq. (8), e = 0) and for a strong circumferential variation in surface 
temperature T = T^ + - T g ) cos 0 at r = R Q (eq. (8)). Consider first the constant 

temperature solution (e = 0) in detail, then the 0-dependent solution. 


Constant Surface Temperature, T = T,j 


The nondimensional form of the energy equation (eq. (6)) becomes 


JL k 2 = 0 
a? V a?/ 


where 


0 = 


T - T ■ 


T - T. 
s d 


which can be solved directly. From appendix F the solution is 


(B 18 ) 


(B19) 
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9= iK) 


(B20) 


The temperature gradient at each surface becomes 


dQ _ (£_ 
A 


(B21) 


d0 1 

£=1+A <? A 


(B22) 


Interface energy balance . - The velocity and temperature distributions have been 
expressed up to this point in terms of an unknown parameter, the dimensionless vapor - 
gap thickness A. The interface energy balance (eq. (13)) is now used to determine the 
value of this parameter. 

First, the radial velocity at the interface must be determined. Substituting the 
values of the constants Cj to (eqs. (BIO) to (B12), and (B16)) into equation (B3) and 
evaluating the velocity at the interface (£ = 1 + A) gives 


/15Wj 
A \4irw* 


1 + A + — A 2 

15 


A cos G 


cp^l + (p + <p^ + 6 (<p^ + <p^)j ^ 


(B23) 


For very small values of A, v^ may be approximated by 


V A«1 “ ' 


W d\ A 2 cose 


Anw*J g(0*) 


(B24) 


where 


w 


* 




s. 


(B25) 


Substituting the expression for v^ (eq. (B23)) into the interface energy balance (the 
nondimens ional form of eq. (13)), imposing the constant-temperature boundary con- 
straints, and solving for A give 


A = f(0) 


(B26) 


This gives rise to a serious problem in the mathematical operation at this time, since 
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A is not independent of 9. Consequently, the similarity solution for this transform 
does not exist in terms of the variable i }/(^, 9) for the constant-temperature case, since 
A had been assumed a constant in the solution of the momentum equations. However, 
the constant-temperature case is of great interest, and the concept of similarity condi- 
tioning is now introduced to attain a solution. This problem does not occur in the second 
case where surface temperature varies as cos 9, as is discussed later in this appendix. 

Similarity conditioning . - From a practical engineering point of view, the answer to 
the similarity dilemma posed herein was to reduce the severity of the auxiliary constants 
given by equation (13). As in many analyses, simplifications are made to the physics 
and then checked for validity. Thus, instead of requiring the interface energy balance to 
hold at every point on the boundary, as required in equation (13), an integral boundary 
constraint of the following form is used in place of equation (13): 


-P XV r (R o + 6) 



r - R o +6 


(B27) 


where 


V r (R 0 ♦ 5) = 


C f [V(R 0 + 6)](R q + 5) 2 sin 9 d0 d$ 


/*277 y-e* 2 

/ / (R. + 6/ sin 9 d<9 d$ 

* / 0 * / 0 o 


(B28) 


Here, V r represents an average radial velocity in the vapor gap. This integral con- 
straint conditions the problem sucn tnat a similarity transformation exists. 

Similarity solution by integral boundary constraint: Nondimensionalizing equa- 
tion (B27) results in the averaged energy balance 


= s — 
A dC 


C=l+A 


(B29) 


where 


S 


Pr 


C p (T c 


t s v 


(B30) 


At this point, 
independent of 9. 


A is an unknown function of 9. 
In seeking an average value for 


If v^ were constant, 
v^ independent of 9, 


A would be 
A, while still 
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an unknown, would also be independent of 9. Therefore, assuming A as some average 
value and substituting the expression for v A (eq. (B23)) into the nondimens ional forms 


v A becomes 



( Wi+ S a2 ) (1 

A y 8?rw* J 

15<p(l +3A + Ha 2 + 2 A' 
\ 3 


(1 + cos 6*) 


| A 4 jg(0*) 


(B31) 


Substituting equation (B31) into equation (B29) along with the expression for the gradient 
(eq. (B22)) and solving for A gives 


a 4 _ / 87rw*S i 
\ W, j 


(l + 3A + — A 2 + 

2A 3 + — A 4 \g(0*) 

\ 3 

5 / 

U + A + i- A 1 

2 Vl + cos 0 *) 

\ 15 

j J 


For A « 1 and G(6*) = g(0*)/(l + cos 9*), equation (B32) becomes 


(B32) 


As r 8 ,m*SG(nl 1/4 (BM) 

L w d J 

Thus, A is independent of £ and 9, but dependent on 9*, and the similarity solution 
exists. 

Approximate solution for small angular displacements: As discussed previously, 
some of the natural physics of the problem were relaxed in order to preserve similarity 
and obtain a separable solution to the differential system. 

In this section, it is shown that, given a value of angular position 0, similarity will 
be approximately preserved for small angular displacements. This requires that to 
C 4 be independent of angular displacement. The constant C 2 is already an absolute 
constant, as is now shown. Substituting the value of C 4 (eq. (B16)) into the expression 
for C 2 , equation (Bll) gives 


C 


2 



(B34) 


Thus, the constants Cj, Cg, and C 4 need only be examined in detail to determine their 
dependence on £ and 9. The variation of Fj (eq. (B2)) can then be determined. 
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As 9 approaches rr/2 (which implies that 6* — n/2), A^ becomes very much greater 
than 1. 

Since C^, Cg, and C 4 are functions of A (see eqs. (BIO) to (B13), and (B16)), 
they turn out to be pseudoconstants which are constrained by sec 9. However, the solu- 
tion to the differential system was based on a separable solution. The constants Cj, 

Cg, and C 4 being functions of 9 clearly violate this premise. Nevertheless, approxi- 
mate similarity would be preserved, if for small angular displacements the constants 
were indeed invariant. The variation with angular position in the pseudoconstants, the 
gap thickness A^, the interface velocity v^, and the pseudofunction Fj is presented 
in figure 5, using the approximate solution for A. 
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Figure 5. -Variation of gap thickness A, interface velocity v*, and pseudo-steady- 
state parameters with angular position (0* - tt/2, £ = <p, ana Q= 3.55X10 5 ). 


As seen in figure 5, for 9 < 30° the pseudoconstants Cj, Cg, and are very 
nearly constant; however, A^, v^, and F vary 3, 4, and 10 percent, respectively. 

For 6 < 86°, they vary as follows: 

Cj < 4 percent; Cg < 3 percent; C 4 < 8 percent (B39) 

However, as evident from figure 5, F changes by a factor of 6. 65, which is clearly a 
violation of the separability hypothesis. Nevertheless, the pseudo- solution points out 
that as A^ increases, v^ decreases; this means that the heat transfer at the interface 
is decreasing, which is in agreement with reality. These errors at small angles are quite 
acceptable considering the physical intuition gained by this approximate technique; how- 
ever, at large angles, the error in F is considerable, and the results are highly ques- 
tionable. 
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The values of and Ag in equations (B33) and (B38) differ only by a factor 


2 cos 9 X 1 / 4 
1 + cos Q*/ 


(B40) 


However, at 6 equal to 60° and assuming 0* - 77/2, both the similarity and approximate 
solutions for A agree. This implies that both solutions can be made to agree in an av- 
erage sense, which, in turn, implies the surface-averaged heat-transfer coefficients are 
about the same 


h A « h g (B41) 

Heat-transfer coefficient . - The heat-transfer coefficient can now be constructed 
from the momentum and energy solutions. Because the approximate and similarity solu- 
tions presented in the last section differ only by a constant when average values are con- 
sidered, only the similarity case is considered here. 

Sphere surface: The heat-transfer coefficient to the sphere is defined as 


h d< T d- T s> = - k f 


r=R, 


(B42) 


Utilizing the nondimens ional forms, h^ can be written as 


k 


o 



Substituting equation (B21) into equation (B43) gives 


h d R o 



(B43) 


(B44) 


Substituting the expression for A given by equation (B33) into equation (B44) gives 
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Nusselt number Heat con- Concentric-surface 
for sphere duction Leidenfrost 


from phenomenon 

sphere to contribution 

stagnant 
fluid 


(B45) 


In dimensional form, the gap thickness can be expressed as 


87ri'k(T d - T s )G(e*)] 1//4 

P dV x 


This was accomplished by substituting the dimensional forms of S (eq. (B30)) and w* 
(eq. (B25)) into equation (B33). In a similar manner, the heat-transfer coefficient may 
be written 



3 T 

Pd v dg** k 

87rR^(T d - T g )G(0*) 


(B47) 


where X has been replaced by X* to compensate for convection effects, see appen- 
dix G. But 


therefore, 




(B48) 


k 3 P d Pgx’ 


R o^ T d- T s) G ( 0 *) 


This equation indicates that the stagnant conduction term k/R Q will begin to dominate 
whenever the radius of the sphere is such that 
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-5- <* 0. 64 
R„ 


P d PS x> 


R 0 P( T d - T s )G(0*) 


1/4 


(B50) 


or 


R Q < 1.81 


M T d - T s )G(0*) 

P d Pg>* 


,1/3 


(B51) 


That is, when the gap thickness becomes larger than the sphere radius, stagnant con- 
duction would dominate; and, of course, the entire solution becomes questionable. 

If at this point the similarity solution is compared to the approximate solution (refer 
back to section Interface energy balance ), the approximate solution indicates that as 0 
approaches v/2 (which implies 0* — ir/2). 


h . = — for 0 — 7t/2 (B52) 

R o 

That is, the heat-transfer coefficient returns to the stagnant flow contribution of a sphere 
in an infinite media. 

Liquid interface: The heat transfer at the interface between the vapor and the sup- 
porting liquid is given by 


h i (T d " T s) = " k — 
1 a s 9r 


(B53) 


r- R o+6 


Nondimensionalizing this equation and substituting equation (B22) into the result yield 


h i R o 


^ - 1 


A(1 + A) A 


(B54) 


Substituting the value for A given by equation (B33) into equation (B47) gives 


h.R„ 

Nu. = -J-H = 
1 k 


w a 11/4 


L8ttw*SG(0*)J 


- 1 


(B55) 
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In this case, in comparison to equation (B45), the Nusselt number equals the concentric- 
surface Leidenfrost phenomenon contribution less that of the axisymmetric stagnation 
flow. 

The approximation solution (refer back to section Interface energy balance ), how- 
ever, shows that as 9* approaches n/2, A becomes large, and the approximation 
used in equation (B54) is no longer valid. The value of is zero at this point. Con- 
ceptually, at 6* equal to jr/2, the sphere is emitting energy to an infinite stagnant 
region; that is, 


h, =A for 9*~- (B56) 

d R 0 2 


while the energy never reaches the supporting interface at infinity; that is, 


h. = 0 for 9* - - (B57) 

This implies that there can be no evaporation at the liquid interface at 0 = n/2 because 
of the presence of the drop. 


Strong Surface-Temperature Variation 

The second problem is that of a surface temperature variation of the form 

T = T d + e(T d - T g )cos 6 at r = R Q (8) 


The energy equation for a strong surface-temperature variation with angular posi- 
tion is solved in appendix F. The solution is 


where 


ae, 

a? 


cos 9 


(1 + 2cp J ) 


c=l W - l)cos 0 ref 


a©„ 




3ej cos 9 


£=1+ A (cp* - l)cos 9. 


ref 


(B58) 


(B59) 
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and is an arbitrary positive constant (i. e. , > 0). 

The temperature gradient at each surface becomes 


0 e -l— L - cos -— for 9* < - 

" <p 3 - l\ C 2 / C ° S0 ref 2 


(F12) 


where 0 C is defined 


Interface energy balance . - The velocity and temperature distributions are ex- 
pressed in terms of the dimensionless vapor-gap thickness A, and again the interface 
energy balance is used to determine A. 

The radial velocity at the interface (eq. (B23)) and the temperature gradient at the 
interface (eq. (B59)) are substituted into the interface energy balance equation 
(eq. (B35)) 


T - T 


®c = 


( T d) 


- T, 


ref 


(F10) 


-v A = S 


d0 

d? 


(B35) 


?=1+A 



(l + A + ±A 2 ^cosJ 


15 


/Lg(e*) 


] 


3Sej cos 9 


1 + 3A + — A 2 + 2 A 3 + 


| A 4 ^ ( (p 3 - l)cos 0 ref 


(B60) 


where 


g(9*) = 1 - 3 cos 9* + - cos 3 0* 
2 2 


(D26) 


Solving equation (B60) for A gives 


a ( 4irw*Se 1 x 

A 4 = ( — W 


1 + 3A + — A 2 + 2 A 3 + 2 A 4 


V W d cos ^ref / 




g(9*) 


(B61) 
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Now for A « 1, equation (B61) reduces to 


A 


4 


47rw*Sejg(0*) 


W d cos 0 ref 


(B62) 


The similarity of equations (B62), (B33), and (B38) is striking; for example, if — 1 
0 re j — 7i/3 and 6 — 77 / 3 , where 6* — ti/2, they are the same. 

Heat-transfer coefficient . - As in the constant- surface-temperature case, the heat- 
transfer coefficient can be constructed from the momentum and energy equations. 

Sphere surface: The heat-transfer coefficient to the sphere is defined as 


V T d 


T s ) = -k 


3T 

3r 


r=R 


or writing equation (B63) in nondimensional terms, 



(B63) 


(B64) 


Evaluating equations (F12) and (B58) at £ = 1 and substituting into equation (B64) gives 


k 


o _ 1 + 2 £_ 

<P 3 ~ 1 


(B65) 


Assuming A « 1 and substituting for A (eq. (63)) in the result, equation (B65) be- 
comes 


h^R 

k 


0 ^ 



1 + 


W d cos 0 ref 


477W*Se 1 g(0*) 


1/4 


(B66) 


H 0 re £ = 60° and — 1, where 6* — n/2, equation (B66) gives the same Nusselt num- 
ber as for the constant-surface-temperature solution (eq. (B45)). 

Liquid interface: The heat-transfer coefficient at the interface is given by 
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Substituting equations (B59) and (F12) into equation (B67) gives 


h.R o i 

- 1 ° = — i — » - - 1 (B68) 

k 3 i A 
<P - 1 


Equations (B68) and (B54) are identical. 

The equivalence of equation (B66) to (B45) and of equation (B68) to (B54) indicate 
that, insofar as heat transfer is concerned, the assumption of a constant-temperature 
surface with a constant vapor-gap thickness is compatible with a strong 0-temperature 
variation with variable gap thickness. 


NONDIMENSIONAL FORMS 


It is most instructive to investigate the heat-transfer- coefficient expressions for 
nondimensional parameters. The nondimens ional forms identify the major, or first- 
order, terms influencing the solution. In equations (B45) and (B66), the term W^/w*S 
appears, which can be expressed as 


w*S 




q 

R oP 2 «c 


Pr 


Cp AT 



p(pl - p)R^g 


2 

M 




(B69) 


The latter part of equation (B69) is the familiar Grashoff number - Prandtl number 
product (or Rayleigh number modified by the ratio of latent to sensible heats), and X is 
replaced by X* 


Ra’ s Ra = GrPr — 

CL AT C AT 

P P 


(B70) 


Substituting equations (B69) and (B70) into equations (B45) and (B66) gives 


1/4 


(B45a) 


Nu d 


Ra' 

p d \ 

60(0*) 

\ p l ' p ). 


Nu d 


1 + 


Ra' 

p<i \ 

cos 6 r ef" 

3g(0T 

\ p l ' p ) 

e l 


-11/4 


(B66a) 


Equations (B45a) and (B66a) indicate that Nusselt number is governed by the usual film- 
boiling relation (Ra*)*^, the interrelation between density differences, and 9*. 
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APPENDIX C 


SOLUTION OF HEAT -TRANSFER COEFFICIENT FOR SECOND 
SIMILARITY TRANSFORM T 2 (WHERE 0*<7T) 


MOMENTUM EQUATION 

The second similarity transform 

0) = f 2 (0(i - cos 6 ) (30) 

is used to convert equation (27) into an ordinary differential equation. This transforma- 
tion has also been used for solving for flow around a sphere. Substituting equation (30) 
into equation (27) yields 


f |(0 = 0 

This is a linear homogeneous equation of fourth order which has a solution 

f 2 (iq=c 1+ c 2 c+ c 3 ? 2 + c 4 ? 3 

Therefore, 



(Cl) 


(C2) 


(C3) 


Fo (C)tan i 

v. - - * (C4) 

u ^ 

The constants in the previous equation must be determined from the boundary con- 
ditions (eqs. (8) and (9)). In terms of the dimensionless variables, these boundary con- 
ditions become 


v 0 (l, 9) = 0 


(C5) 

(C6) 
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v 0 (l + A, 9) = 0 


Vj» (1 > 9) = 0 


(C7) 


v^(l + A, 0) = v A 


(C8) 


where 


A = -L (C9) 

R o 

Applying these conditions to equations (C3) and (C4) allows the four constants to be 
expressed in terms of one constant of integration 


C-i = — (1 - 3 <p) 
1 2 

(CIO) 

^2 = *^4^ 

(Cll) 

3C 4 


c 3 = - (1 + (p) 

(C12) 


2 


where 


cp = 1 + A (Cl 3) 

The dimensionless gap thickness A is assumed at this time to be an unknown con- 
stant. However, as shown in appendix B (eq. (B27)), this requires that the interface 
energy balance be satisfied on an average over the heating surface rather than at every 
point. 

The constant C 4 can be found by satisfying the static force balance constraint. 
First, however, the pressure distribution in the vapor gap must be found. Substituting 
equations (C3) and C4) into equations (20) and (21) and solving for the pressure distribu- 
tion give 

Pl^i Pi = "X-L+JP) - 6 In £ + 12 In /sec + — (C14) 

C 4 s V 2 ) C 4 

The constants C 4 and are determined from the static force balance in appen- 

dix E. This gives 
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(C15) 


W 


C 4 = “ 


S^w 51 


cos 2 0* - 2 cos 0* + 1 


3 r P 


— = + 3(1 + cp)- 12 In (sec — 


C 4 C 4 


V 2 


(Cl 6) 


Therefore, utilizing Cj, C 2 , C 3 , and C 4 , the velocity and pressure distributions 
are known functions of cp, where cp = 1 + A. 

Next, the energy equations are considered. Afterwards, the solution to the energy 
and momentum equations are combined in the interface energy balance to obtain a solu- 
tion for the heat-transfer coefficient. 


ENERGY EQUATION 


The energy equation is solved in appendix F, and the general solution for a cosine 
variation in surface temperature of an arbitrary amplitude is (eq. (F5)) 


e = 2 (l - + — S — (r . £?\ 

aV v „ 3 -i\ t 2 / 


COS 


(C17) 


where 


T - T ■ 


0 = 


T s - T d 


(CIS) 


The temperature gradient at each surface becomes 


30 

ac 


C=i 


= ^ (1 + e cos 0) 
A 


(C19) 


a© 

a? 


1 |i e cos 0 

?=* A n i + a? 


(C20) 
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which for small values of e reduce to equations (B21) and (B22), respectively. Also 
note that where e cos 9 » 1, which corresponds to the strong- variation case and 
6 * < 7t/2, equations (C19) and (C20) are similar to equations (B58) and (B59), respec- 
tively. 


INTERFACE ENERGY BALANCE 


The velocity and temperature distributions are expressed in terms of two unknowns, 
A and 9*. The interface energy balance is used to determine A in terms of 9*. 

Substituting the values of constants to Cg into equation (C3) and evaluating v^ 
at q> give the interface velocity 

+ 3 - <pj (C21) 

where is a function of 0* and the parameter W^/w* (see eq. (C15)). 

Nondimens ionalizing equation (13), the interface energy balance becomes 


= _ 1 ( 1 - 3<y 


<P 


v 


A 



(C22) 


where 


1 

r c P (T d - v 

Pr 

x* 


(B29) 


and X has been replaced by X* to accommodate some effects of convection (see ap- 
pendix G). 

Substituting equations (C15) and (C21) into (C22) gives, for <p & 1, 


A 


4 


( V* S V cosV - 2 cos 9 * + 1) 

V w d / 




(C23) 


The problem here is, as in appendix B, the value of 8* and how to assign e cos 9 
so as to maintain similarity. The latter is solved by assigning 
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€ = 0 


(C24) 


The value of 6 * is determined from the balance of forces at the interface (see appen- 
dix H). 

HEAT-TRANSFER COEFFICIENT 
Sphere Surface 

The Nusselt number for the sphere is given by 


h d R o _ 00 

k 0? 


(B43) 


s=i 


Substituting equations (C19), (C23), and (C24) into equation (B43) gives (assuming 
A « 1) 


h d R o 


= 1 + 


W. 


[6j7W*S(cos^<9* - 2 cos 0 * + 1)J 


1/4 


(C25) 


where 9 * is determined in appendix H. 


Liquid Interface 


The nondimensional form of equation (B53) gives the Nusselt number at the inter- 


face 


h i R o _ dQ 

k ac 


(C26) 


S=cp 


Substituting equations (C20), (C23), and (C24) into (C26) gives 


k A (p A 


(C27) 
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h.R^ 

-L-° = -1 + 


W. 


|_3rw*S(cos^0* - 2 cos 0* + 1)J 


1/4 


where, again, 0* is determined in appendix H. 


(C28) 


NON DIMENSIONAL FORMS 

In appendix B the basic nondimens ional terms governing the solution were identified; 
these parameters and W^/w*S are given by equations (B69) and (B70) as 


W 


w*S 


d - 4 >Ra'_f! 


Pj - P 


(C29) 


Substituting equation (C29) into equations (C28) and (C25) indicates that heat transfer is 
influenced by the usual Ra’*^ and the interrelation between the density ratio 
Pd /(pj - p) and 9*, as was shown in appendix B. 


d o 


= 1 + 


- Ra' ( Pd 


\Pl ~ P ^ 


1/4 


(C25a) 


h.R^ 

-L_2= -l + 


— Ra’ I 


p d 


" P ^ 


1/4 


(C28a) 


where 


f(0*) = cos* $* - 2 cos 0* + 1 


(C30) 
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APPENDIX D 


STATIC FORCE BALANCE FOR TRANSFORMATION Tj (WHERE 0*< tt/2) 

The pressure distribution beneath the spherical drop comes directly from the solu- 
tion of the momentum equations and is given by equation (B15) (repeated here for con- 
venience) as 


/ 2C 2\ 3 

p(C, e) = -cos 0^2OC 4 ? + -yj + X 


(B15) 


where C A and C 9 are related by equation (Bll). The purpose of this appendix is the 

~z Ci O 

determination of the two unknown constants C 4 and °C by the application of boundary 
condition equation (11) and the static force balance given by equation (12). 


PRESSURE DISTRIBUTION IN VAPOR GAP 


Evaluation of equation (B15) at the boundary defined by equation (11) gives 


p(l, 9*) = p a = -cos 0*(2OC 4 + 2C 2 ) + X 


(Dl) 


where the dimensionless atmospheric pressure p is given as 




pu 


*2 


(D2) 


Solving for the constant C gives 

3 C = p a + cos 0*(2OC 4 + 2C 2 ) (D3) 

Substituting equation (D3) back into equation (B15) gives for the pressure distribution 


p(C, 6) = 20C 4 (-? cos 6 + cos 9*) + 2C 2 | 


-cos 


+ cos 9* 1 + p r 


(D4) 


The pressure and velocity distributions which relate to the shearing stresses are 
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now substituted into the static force balance (eq. (12)) to determine the remaining con- 
stant. The static force balance is divided into four integrals so that the contribution of 
each term can be more easily visualized. Figure 2 shows a schematic of the pressure 
and shear forces acting on the sphere. 


FORCES ACTING ON SPHERE 


Pressure Force 


The component of pressure force within the gap given by equation (12) is rewritten 


as 


/ 2n j»8* 

/ P(r, 8) | R cos 6 R* sin d d 9 d$ 
u y 0 o 


, 2 


(D5) 


To nondimensionalize the weight as well as the components of force, the parameter w* 
is defined as 


pu* 2 R/? 
w* 2 


(B25) 


Dividing equation (D5) by equation (B25) gives 



However, since 


P(r, 9) _ P(C, 9) 
w* d 2 


(D6) 


m 


equation (D6) becomes 


I 


1 



cos 9 sin 9 d 9 d$ 


(D8) 


47 



Evaluating this integral gives 


I 


1 


4?r £(10C 4 + C 2 )^- i i cos 9 * 


- cos 3 0*\ + — (1 - cos 2 0*) 
6/4 


Tangential Shear Force 


The component of force due to tangential shear in the static force balance 
as 



T re sin 0 lr=R R o sin0d0d * 


Nondimens ionalizing this equation gives 


h = — = 


w n 



- sin el n sin 0 d0 d<f> 
w* o 


But 


T t > e ~ 


're 


w g 


c r 0 w*L ac\c / ? 


However, 


* 

MU _ M 

g c R 0 W * PU*R 0 rJ 


1 _ 1 

ReR n R n 

O O 


since the Reynolds number is identically 1. Thus, 



sin 2 0 d0 d$ 

?=1 


(D9) 


is given 


(DIO) 


(Dll) 


(D12) 


(D13) 


(D14) 


48 



Substituting the values of v 0 (l, 9) and v ? (l, 0) into equation (D14) and integrating yield 




+ CA 


cos e 


(sin 2 0* + 2) 


-f] 


(D15) 


Normal Shear Force 


The normal component of shearing force acting on the sphere is given by 

•2 n /'O* 


pin p 

•'O *'0 


r rr cos 0| r _ R R q sin 9 d 9 d$ 


(D16) 


Nondimens ionalizing this equation gives 



But 


(D17) 


r - T rr _ -2juu* _ -2 
^ w* w*g c R Q a? r 2 ac 


(D18) 


Therefore, equation (D17) can now be rewritten as 

■»2n p9* 


I 3 = 


pin pb 
*'0 •'0 


0V^ 


cos 6 sin 9 d0 d$ 


(D19) 




and since dv^/dt; |^ =1 = 0, 


V 0 


(D20) 


Ambient Pressure Force 

The ambient pressure acting on the upper portion of the sphere becomes 
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(D21) 


rl-n rQ* 

^-1 / P o 


'0 7T 

Nondimensionalizing and integrating gives 


cos 0R* sin 6 d0 d$ 


I 4 = 77P a (l - COB 2 6*) 


SUMMATION OF FORCES 

Summing these dimensionless forces acting on the sphere gives 


W, 


w 


_ h + h ' h " t 4 


or, after substituting in the values of to 1^ 


W, 


Anw*C / 


1 

(.0*^ 

1-4 +-j] 

(l - cos 

8* + - cos 3 9*) 

3 

V C 4/ 

1 \ C J 

V 2 

2 / 


Now Cj and C 2 are known values of C 4 ; however, 0* is still an unknown. 


W, 


C 4 = 


. 87TW* 


1 + <p + 4a/ 


1 + cp + <p 2 + 6(cp 3 + (p 4 )]g(0*) 


where 


g(0*) = 1 - - cos 9* + - cos 3 0* 

2 2 

Based on experimental observations of water drops floating on liquid nitrogen, i 
that 6* is near v/2. Thus, for 8* near 77 / 2 , equation (D2 5) simplifies to 

1 + <£ + 

1 + cp + <p 2 + 6 {cp 3 + cp 



(D22) 


(D23) 


(D24) 


(D25) 


(D26) 

seems 


(D27) 
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APPENDIX E 


STATIC PRESSURE BALANCE FOR TRANSFORMATION T 2 (WHERE 0*<tt) 


PRESSURE DISTRIBUTION IN VAPOR GAP 

The pressure distribution within the gap is obtained from the solution of the momen- 
tum equations (appendix D). It is given by equation (C14) as 


p(^_0) = -3(1 +J£) _ 6 l n £ + 12 In (sec -'j + — 
C„ C \ 2/ C„ 


(C14) 


where 6 C and are constants. The purpose of this appendix is to determine these 
constants subject to the boundary conditions 


Weight of sphere ~ f Pn • 
* / A 


dA 


(El) 


and 


(C,0) - (i,e*) 

P-P a 

Applying condition (E2) to equation (C14) gives 


(E2) 


— = -3(1 + cp) + 12 In (sec — \ 

C„ V 2 / 


2 / + C, 


(E3) 


or 


+ 3 (i + (p) - 12 In (sec \ 
C, C, \ 2 / 


c 4 c 4 


(E4) 


Substituting equation (E4) into equation (C14) gives 


P(C, 0) - P a 


- -3(1 + <p)(l - lj - 6 In C + 12 


In 


' sec — 


isec -2- , 


(E5) 
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The remaining constant may now be determined from the static force balance 
(eq. (El) or eq. (12)). 


FORCES ACTING ON SPHERE 


Pressure Force 


The force component due to pressure is given by 

•27 r rO* 


f* * 9 

L = / / P(r, 9) | R cos 9 Rr sin 6 d 9 d$ 

y Q y 0 ° 


In nondimens ional form, equation (E6) becomes 

'277 y.0* 

W’ •'o ^0 


i /2 tt ,9* 

Ij = — = J J p(£, 0)j^ =1 cos 9 sin 0 de d$ 


(E6) 


(E7) 


where w* is given by equation (B25). Evaluating this integral using equation (E5) gives 


L = +277C, 



12 In 


sec - 


, sec 


cos 9 sin 9 d0 


(E8) 


or 


\ = 2t7C 4 


a (cos 2 0* - 1) - ^ cos 2 9* + 3 cos 9* - - 
2 


2C, 


(E9) 


Tangential Shear Force 


The tangential shear component is given by 


277 r9 




0 0 


r r0 sin0| r=R R 2 sin0ded4> 


(E10) 
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In nondimens ional form, equation (E10) becomes 


2 7T z.0* 


A../"/ 


*o - - / / r» a sin 0 L * sin 0 d0 d$ 

2 w * % ^’ d ?=1 


and substituting for equation (Ell) becomes 

*2n /» 0 * - 

k t 


I 2 = 




sin 0 d0 d<f> 


Evaluating the integrand and recalling that 


v 0 = 


Fl, tan — 
_ 2 


where F«(l, 0) = 0, and 

F (1, 0) = 2Cg + | = [-3(1 + cp) + 6]C^ 


then 


l 2 = 67rC 4^ “ 


S"!*! . cos e* + l S 


(Ell) 


(E12) 


(E13) 


(E14) 


(E15) 


Normal Shear Force 

The normal shear force acting on the sphere is given by 

'2tt /* 0 * 


/'4TT 

*3 = / J T r 

•'A » A 


• rr cos e p R sin0d0d$ 
'0 -0 R o 0 


(El 6) 


Now nondimensionalize equation (El 6) using equation (B25) and recall that from equation 
(D18) 
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(E17) 




rr 


w 


R 


3v^ 

2 


Evaluating this derivative at £ = 1, 




= 0 


leads to the trivial result 


l 3 ~ 0 

Ambient Pressure Force 

The ambient pressure acting on the upper portion of the sphere becomes 

/'27T „ 

i A = / / p„ cos sin 6 d0 d$ 

4 y 0 y n a 0 

Nondimensionalizing and integrating 

h = ^a ^ 1 “ cos20 *) 


SUMMATION OF FORCES 


Summing the forces acting on the sphere gives 


^d 

w* 


- h + h ~ 



Substituting in the values of equations (E9), (El 5), (E19), and (E21), and assuming 
be small gives 


(E18) 


(E19) 


(E20) 


(E21) 


(E22) 
A to 
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Now, solving for C^, 


w; 

cos^0* - 2 cos 0* + 1 


Equation (E5) may now be written as 


P(C, #) - P a 



e* 

w d 

-3(1 + <p)Q- - lj - 6 In £ + 12 In 

cos — 
2 

3ttw* 

9 

cos - 

f(0*) 


2 



where 


f(0*) = cos^0* - 2 cos 0* + 1 


(E24) 


(E25) 


(C30) 


Now that is known, the velocity profiles are also known. 


APPENDIX F 


SOLUTION OF ENERGY EQUATION (NO CONVECTION) 


Nondimensionalizing the Laplace form of the energy equation with no $ dependence 
due to symmetry gives 


V 2 0 


_L _L U 2 i®\ , 1 

? 2a ?V ^slne 



(FI) 


where 


0 = 


T s" T d 


(F2) 


SOLUTION OF V 2 0 

The boundary conditions are defined such that the temperature at the evaporating 
interface is at the saturation temperature 

T = T_ (evaporating interface) (9) 

s 

and the temperature of the opposite surface follows a cosine distribution, 

T = - AT cos 0 (nonevaporating surface) (8) 

In nondimensional form, the boundary conditions given by equations (8) and (9) for an 
evaporating supporting liquid become 


e{cp,e) = l 


(F3) 


0(1,0) = -e cos 0 (F4) 

where e represents some fraction of the overall temperature difference (i.e. , 

AT = e(T d - T s )). 

Solving equation (FI) subject to equations (F3) and (F4) gives the temperature distri- 
bution in the vapor gap 


56 


0 = £ (l - + — - — K - £- jcos 0 (general solution) (F5) 

* \ V cp*- 1 V c 3 / 

For sufficiently small e, equation (F5) reduces to the radial conduction solution with the 
constant-temperature boundary condition 



0(1,0) = 0 

(F6) 


0(<M) =1 

(F3) 

q = £ ( l _ I\ 

(radial conduction dominant) 

(F7) 




SOLUTION OF V 2 0 c , 0* <nl2 

The solutions (eqs. (F5) and (F7)) are not dependent at this point on the position of 
the interface angle 8*. However, one important case to consider is a solution of equa- 
tion (FI) for 6* < tt/ 2 subject to a different but similar set of boundary conditions 

e(l,fl) = € co -- (F8) 

COS 0 ref 

O c (cp,8)=0 (F9) 

where e is an arbitrary positive constant, e >0, which is at the disposal of the experi- 
mentalist and/or theoretician, and 


0 


c 


T " T g 



(F10) 


0 < 0 < 8 * <- (Fll) 

ref 2 

(Note: 0 c * 0. ) Solving equation (Fl) subject to the boundary conditions (F8) to (Fll) 
gives the temperature distribution in the vapor gap 
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€ 


cos 9 


for 9 * <- 

2 


(F12) 


0 c = ' 




[S - 2 - 


? 2 / cos 0 ref 


Equation (F8) represents a 9 -dependent variation in surface temperature confined to 
those cases where the floating interface angle 9 * < ir/2. Furthermore, a form similar 
to equation (F12) can be obtained from equation (F5) by assuming e sufficiently large, 
which reflects the strong 0 -dependent surface-temperature condition. 

The advantage of equation (F12) is, as pointed out in appendix B, that for the special 
case of 0* < jt/2 a similarity solution is attained. One disadvantage of permitting a 
0-dependent surface -temperature variation (eqs. (F4) and (F8)) is that neither e nor 
0 ref is known a priori and must be assumed or dictated by the experiment. The major 
advantage of the constant-temperature solution (eq. (F7)) is its simplicity, but it requires 
the assumption of the average or representative surface temperature. 
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APPENDIX G 


EFFECT OF CONVECTION ON ENERGY EQUATION 

While the convection effects are analyzed for the case where i = Fj sin 0 and 
6* < -n/2, the results are not expected to significantly deviate for the case where 

^2 = F 2 ^ “ cos 

If the flow in the 0 -coordinate is assumed to approximate the fully developed case 
(i. e. , 9T/ 90 is negligible and radial conduction dominates) and viscous dissipation to be 
small, the energy equation reduces to 


V ± ± h 2 *e\ 

5 3? Pr ,2 3C \ 3?; 


where 


T ' T d 

T s" T d 


If the convection terms on the right side of equation (Gl) were neglected, the energy 
equation would reduce to the pure conduction equation given by equation (F7) in appendix F 
of this report. 

Equation (Gl) may be rearranged to give 


v ? Pr _ 1 iL 
^2 ^2 acL V ^ 


Integrating equation (G3) over $ , 0 , £ gives 


/» /»2n /»0* 


$=0 «/ 0 = 0 


v c Pr\ 2 

-i — \ r sin 0 d0 d$ d? 


»2jt /»0* 


ln(C 2 — )lc 2 sin 0 d0 d$ d£ (G4) 
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Substituting for Vj» using equation (B3) and performing the integration indicated by equa- 
tion (G4) gives 


-Pr(cos 8 * + 1) 


/ 


d? = In /*C ? 2 


(G5) 


which represents the variation in 0 due to ? where 0 has been averaged over 8 and 
<$. Solving for 30/3? gives 


30 

3 ? 



exp 


-Pr(cos 8 * + 1) 



Integrating equation (G6) gives 



(G6) 


(G7) 


The boundary conditions of temperature (eqs . (8) and (9)) can be rewritten in dimen- 
sionless form, where e — 0 and T — T^, as 

0(1,0) = 0 (F6) 

e(cp,8) = 1 (F3) 

Applying these boundary conditions to equation (G7) gives 



Therefore, the temperature distribution becomes 
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exp 


-Pr(cos 0* 


+1 , fm 

J\ y 


dy 


0 = 


/ 

T 2 

- U i 

t/l 



S*<P 

/ exp 

-Pr(cos 0* + 1) f Zixl dy 

A 


1 — 

l y _ 

dr 


(G10) 


However, expanding the exponent in equation (G10) in a power series about 1 and per- 
forming the integration give 


0(0 



Conduction Convection 


correction 


(Gil) 


where 


A(0A) 



(G12) 


The first term in equation (Gil) is identical to the solution for pure conduction, 
while the A(£,a) term accounts for convection effects. Conduction is assumed to dom- 
inate; thus A(£, a) is taken to be 1. However, the range of applicability of this assump- 
tion is now considered. 

Dividing equation (G10) by equation (Gil) gives 
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The parameter A was evaluated numerically (using double -precision arithmetic) by 
employing Simpson's rule to evaluate the integral in equation (G13). As a specific exam- 
pie, a plot of the numerical results is shown in figure 6 for W d /8jrw* equal to 0.472x10 
(water drop on liquid nitrogen). The temperature distribution is shown to be greater than 
predicted from the conduction solution. This is expected, since the specific heat of the 
vapor prevents it from dropping to the value predicted by conduction alone. As shown in 
this figure, the conduction solution is within 10 percent of the actual temperature profile. 

The slopes at the sphere surface (C - 1 = 0), and at the supporting liquid surface 
(C - 1 = 0.04), are a measure of the heat-transfer coefficient since the heat-transfer 
rate is, of course, proportional to the temperature gradient. As shown in figure 6, the 
gradient at the sphere surface is greater than that at the supporting liquid surface; that 
is, the heat-transfer coefficient from the sphere is greater than that to the supporting sur- 
face. This is expected: since all the heat leaving the sphere does not reach the support- 
ing surface, some of it is convected away in the vapor stream. 

The analytical expression for the temperature gradient at the drop surface can be 
found by differentiating equation (Gil) and evaluating the derivative at 5 = 1. This gives, 


9C 


= (1+A- 1 ) 


C=1 



Conduction Convection correction term 


(G14) 


where 


SUM 


^ ^ j-Pr(cos 6 * + 1) 

fel 



(G15) 
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The magnitude of the convection term was evaluated numerically. This term required a 
10-percent correction in the temperature gradient for the case where Pr = 0. 7 and 
A = 0. 04. However, at the expense of accuracy and for increased simplification, the 
closed-form conduction solution was used in the analysis with A replaced by A*, 


A* = 






(G16) 


which tends to correct for the convection omission. 
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APPENDIX H 


SHAPE OF LIQUID-VAPOR INTERFACE 

It has been shown that the sphere is supported by shear and pressure forces acting 
at the vapor-sphere interface, see appendixes D and E. In this appendix, it is demon- 
strated that the fluid possesses sufficient surface strength and buoyancy to restrain the 
shear and pressure forces acting at the vapor-liquid interface. 


FORCES AT INTERFACE 

Within the supporting liquid, there are two forces to be considered, buoyancy and 
surface tension. Buoyancy is composed of two terms: (1) fluid displacement by the 
sphere to an angle of 6 * and (2) fluid head as a result of fluid displacement by surface 
curvature. Each of these forces is now considered separately. 


Surface-Tension Support 

If a load is placed on an axisymmetric shell or membrane structure (see fig. 7), it 



Figure 7. - Loading of axisymmetric shell or membrane. 
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is supported at the circular interface, or load circle. Here the effective load is supported 
by a component of surface tension 

F a = 27rR o a sin 2 e* (HI) 

F a may also be found by simply integrating the lift component of the surface-tension 
pressure rise over the shell, that is, 

cos 6 sin 9 dd d$ (Hla) 

by Sphere 

The "pressure forces" acting on the spherical interface are illustrated in figure 8. 




Figure 8. - Forces acting at vapor-liquid interface. Fluid displacement by sphere. 


The buoyancy resulting from these forces is obtained by integration; however, it must be 
noted that 9* and Z Q , while constants, are also unknown and must be determined from 
the force balance and interface curvature 

/ 2v r9* r2n 

J 0 PjeosfldA J o J o P v cosf>dA 2 (H2) 
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F b -2» 


I \ h l p ~ + [ Z o + V cos 0 " cos e *)]Pj — 1 cos 0 R q sin 0 d0 

Jn [ g C g cJ 


- 2tt 


pv-e* 

I <P — [ h j + z 0 - R 0 ( 1 + cos 0*)] + P — H Q (1 - cos w)l sin w do> (H3) 

J 0 l S C g C i 


Integrating equation (H3) gives the buoyancy in terms of 0* and Z Q as 


F b = 7T(p z -p)JL R 2J 2 R o 


p l + p COS 0* cos^ 0* 


l_3(P i - P) 2 


6 


+ Z Q (1 - COS^ 0*)J| 


(H4) 


If (p^ - $ — Pj , which is true for most film boiling cases, [(p^ + p)/(p^ - pj] — 1. 


Fluid Displacement by Interface Curvature 


As 0 increases beyond 0*, the vapor gap between the sphere and the liquid becomes 
very large. This implies that the pressure acting on the interface approaches atmos- 


A COS d$ = A“\ 




Figure 9. - Forces acting on a thin shell or membrane. 
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pheric; however, the pressure on the liquid side of the interface is above atmospheric 
and must be balanced by curvature of the interface . 

Consider an element of the axisymmetric surface of figure 8 as shown in figure 9 . 

Let the surface forces per unit length in the 3>- and 6 -directions be designated as 
and Ng, respectively. Let there be an X-Y plane tangent to the surface at some point P. 
Let Z be the coordinate co linear with the radii of curvature and Rg. 

The force A is a result of a change in and with respect to 

A= ( R 3 + ^' W ) de ( N * + ? Ld4 ) (H5) 

Resolving the forces as shown in figure 9 and neglecting second-order terms results in 
the force balances for Y- and Z-forces. 

Y-forces . - The Y-force balance equation is 

2RoN a 

— 2 - 2 . d$ d6 - (RjNg d0 d$) cos (tt - $) + Y p d$ d6 = 0 (H6) 

3 $ 

where Y p is some surface force acting in the Y-direction. 

Z-forces . - The Z-force balance equation is 

-R 3 N $ de d# + (RjNg d$ dd) sin (jt - $>) + (P^ - P^RgRj d$ d9 = 0 (H7) 

where (P^ - P y ) dA is the arbitrary force acting in the Z-direction. 

If Y p s 0 and P - P y = AP, equations (H6) and (H7) become 

0R,N. 

— 2 - 2 . = -R,N„ cos $ (H8) 

3$ A p 



N fl 

- sin $ 


= AP 


(H9) 


Now R 3 = Rg sin ( 7 r - <t>), and substitution into equation (H9) would lead to a description 
of the surface in terms of the radii of curvature, R^ and Rg. 

For this report, we consider a surface of constant strength, = o = Constant, 

being depressed by a nonwetting object; that is, a sphere of water depressing liquid ni- 
trogen in the film-boiling regime. The problem is to solve these equations where Z Q 
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and 6* are unknowns. This requires an iterative -matching solution where the first step 
is to solve equations (H8) and (H9) subject to a given Z Q and 6*. The second step is to 
determine if such a 0* and Z Q will yield sufficient strength to support the sphere. Al- 
though in nature the solution is unique, the computer results may not be unique in that a 
family of approximate solutions may exist near the unique solution. (This could also in- 
dicate that nature is not unique.) Recasting equations (H8) and (H9) by substituting 


u = sin $ 


and nondimensionalizing by the characteristic length 

t 2_ ^ -P)g 


°*c 

which arises in a natural way from equation (H9), give 

Du = - z 

x x 

D z = tan $ = H 

l/i 2 

± f 1 - u 


(H10) 


(Hll) 


(H12) 


(H13) 


where 


z 



(HI 4) 


x 



(HI 5) 


For 6* > n/2, tan $ > 0, and the sign in equation (H13) is positive; for 0* < tt/2, 
tan $ < 0, and the sign in equation (H13) is negative. 

Anticipating the computer solution, the increments are written 


l j + l = u j +(D X U > Ax 

(H16) 

; j + l = z i +(D x Z ) Ax 

(H17) 
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Computer solution. - As the initial conditions are assumed to be known (i.e. , Z Q , 

6 *, and the slope of the interface), a fourth-order Runge-Kutta routine was applied to 
equations (H12) and (H13) subject to equations (H16) and (H17). (The Runge-Kutta routine 
was developed by Paul Swiegert and is obtainable on loan from Lewis. ) The resulting 
axisymmetric profile, shown in figure 10, represents the fluid displaced by interface 
curvature. The problem now becomes one of determining the Z Q and 6 * which not only 
result in a balance of forces but yield sin $ = 0 at Z/L = 0. 




(b) Schematic model for interface configuration. 

Fiyure 10. - Interface configuration for water sphere on liquid nitrogen. Radius of sphere, R 0 = 0.3 centimeter; 
density, - p * O.S gram per cubic centimeter; sphere density, p d = 1.0 gram per cubic centimeter; surface ten- 
sion, o = 8.8 dynes per centimeter. 
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BALANCE OF FORCES 


The forces for a given z Q and 0* must balance the weight of the sphere. This re- 
quires an iterative procedure supplying new z Q and 9* values to achieve the desired 
balance 

1 np s — Rq = F„ + F (H18) 

3 g c 

Herein the solution was attained only for the case, 

R Q = 0. 3 cm 


P z - p = 0.8g/cc 
P s = 1.0g/cc 


a = 8. 8 dynes/ cm 


which is shown as figure 10. While the forces are balanced to less than 0. 3 percent, the 
exact solution was not obtained. The interface is bounded by the given curves: one where 
z — 0 but sin $ > 0, and one where sin $ — 0 but z > 0; and the interface was inter- 

1 ° 

polated from these curves. In this case, 6* is approximately 113-g , and the initial head 
is approximately 0. 148 centimeter. This case was chosen because it is representative 
of the floating spheres of table I, and visual examination of these spheres, along with 
motion pictures, appears to confirm the prediction (see motion-picture supplement 
C-267). 

At the time of publication of this report, a more elegant technique for determining 
an axisymmetric interface was found in the paper by Huh and Scriven (ref. 12). 
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APPENDIX I 


FREEZING TIME: ANALYTICALMODEL 

An experimental check on the theoretical solution requires an estimation of the time 
required to freeze a sphere of liquid. However, there are many facets of the freezing 
problem that we do not understand. For example, the convective heat-transfer coeffi- 
cient from the top of the sphere, as well as internal convection, is unknown. If sub- 
cooling occurs (and up to = -40 K is possible), the solid will not form at the 

freezing temperature, which completely alters a heat balance. The compressive effects 
of freezing on the interface freezing temperature and interface growth rate are not known; 
the effects of cracks, ice "worms,” and the 4° C inversion point are also unknown. To 
what extent does eccentric cooling and freezing alter a concentric-sphere freezing 
analysis? The effect of a pseudo -steady- state analysis, as compared to a complete tran- 
sient analysis, is difficult to assess. To what extent is the metastable Leidenfrost phen- 
omenon influencing the results ? If a sphere followed the metastable Leidenfrost line 
(ref. 4), floating times would be greatly extended. 

Each of these unknowns constitutes a report in itself. As a first-order cut at the 
problem, a pseudo-steady-state model is assumed; the model is illustrated in figure 11. 



Figure 11. - Schematic model of freezing drop. 


Obviously omitted are such effects as subcooling, convection, compression, 4° C inver- 
sion, eccentricity, the transient nature, and cracks. The equation depicting the energy 
distribution within the concentric solid region is 



d T 

solid 9r 


k solid J_ ( r 2 
r 2 3r V dr/ 


(II) 
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(While it must be readily admitted that freezing does not occur in a concentric manner, 
the inclusion of the noncencentric effects is difficult and seems unwarranted.) However, 
if it is assumed that V r is small, energy is transported by conduction, and equation (II) 
can be written 


±/c 2 iIUo 
35 V 35/ 

subject to the boundary conditions 

T - T 

5 = 1, r = — = 0 

Tr ~ 


C = i-/3, r = l 


(12) 


(13) 


where 


0 = 



The solution of equation (12) subject to equation (13) is 



The heat flow out of the sphere may be written as 

\ ’ h d< T o - T s> * h ice< T f - T o> 


Solving equation (16) for T Q gives 


T o = 


T s h d +h ice T f 
h ice +h d 


and using equation (17) in determing yields 


(14) 


(15) 


(16) 


(17) 
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Q r = q r A = 4„R 0 h d 


T f T s 


1 +- 


h. 


ice 


(18) 


From equation (15) 


ar 


?=i 


= £ - i 

fi R r 


- 1 


(19) 


Using equations (19) and (16) results in 


-k. 


h 


ice 


V T o 


R 


ice dT 

3? 


k. 


ice 


\ r f > 


(HO) 


Substituting equation (110) into equation (18) yields the following form for Q R 


Q R “ 


4 * R p( T f - T s> h d 

/R \ R h , 

1 + - A-2J. 


(HI) 


[ f 


k s 


ice 


The interface energy balance gives the advance of the solid surface per unit time, 
which will give the desired freezing time 


~ p ice y 


dr f _ ~ k ice (T f ~ T o> 3 r 
dt R„ 9? 


- V T d " T f> 


( 112 ) 


?=(r f /R 0 ) 


Here T H represents a suitable time-averaged value, that is, 




+ T, 


T d = 


(113) 
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When equations (15), (110), and (111) are substituted into equation (112), the interface ve- 
locity becomes 

r -\ 



By defining the following constants 

1 


C = 1 


Rh, 
o d 

Ir, 

^ice 


2 c - Rphd 

k ice 


(115) 


(116) 



p- vR 
'me' o 


(T„ - T f ) 


4 C _ h d< T J - T s> 
p ice yR o 

equation (114) may be written in a simplified form as 


(117) 


(118) 


/ 


1 - 


1 C 3 Cx 2 + 2 C 3 Cx 


— W/ 

+ 4 cj J * 


5 c dy 


(119) 


where x = r f /R Q . Equation (119) has two solutions, depending on the sign of the charac- 
teristic 



- 4 1 C 3 C 4 C 


Integrating equation (119) gives the freezing time as 


( 120 ) 
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t = 4 


W r f_A + 4 c ln 

2 J C 3 C ^ 

V R o/ 

+ 2 C 3 C - 

2 1 C 3 C + 2 C 3 C + 

3 c\ R o ) 3 C^ 

2 1 C 3 C ( — ^ 

\ R 0 y 

j + 2 C 3 C 

2 ! C 3 C + 2 C 3 C - y^ 


--L(— L-l^d- — 2 C I arc tan ^ 

3 cJT ' 


3 C\ R 0 


2 1 C ( — ) + 2 C 
R 


- arc tan 


3 C(2 *C + 2 C) 


R 


C R >° 


( 121 ) 


for C R < 0 


If the effect of the liquid core is neglected (i.e. , the time required to remove sen- 
sible heat is negligible relative to the latent heat removal time), the solution of equa- 
tion (114) for T d = T f becomes 


, _ <W R o f, 
" 3h a< T d - T s > l 



( 122 ) 


In comparing equations (121) and (122) it would seem well to check the necessity of 
including transport to the liquid core in the interface energy balance (eq. (114)). It is 
also evident that, because of the pseudo -steady -state assumptions, the sensible heat of 
the liquid core could be extracted prior to the freezing time predicted by equation (121). 
The actual time would be between that predicted by equations (121) and (122), depending 
on h^ and T d - T^ in equation (112). 

As the surface temperature cools, the freezing interface advances toward r = 0, 
and the sphere begins to behave as a solid losing heat to the environment. The govern- 
ing equation becomes 


h ,A(T - T, 


■> ’ -Kle' 


0(T 


T s> 


3r 


(123) 


where 


T = t - T, 


(124) 


And at t = 0, the surface temperature T(0) as well as the elapsed 
unknowns. The solution of equation (123) becomes (for /(pCp^ 


time t, are, as yet, 
= Constant) 
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= exp 


(125) 


T - T s 
T(0) - T g 


h d A 



T 


Solving for the time r, 


T 





T(0) - T g 


(126) 


In order to determine T(0), the slopes of equations (121) and (123) are matched and 
the interface location r^/R Q , the time r^, and the temperature T(0) can be determined. 



eq. (121) 




eq. (126) 


(127) 


The change in surface temperature with time is, using the chain rule and letting T = T , 


D t (T - T g ) = 


< T o - T s) 
ice 


D (r f /R 0 ) h ice 


Id. 



(128) 


where from equation (110), 


D (rf/ R 0 ) hice 



and by manipulating equation (17) and substituting equation (16), 


T f - T o 


ice * h d 


(129) 


(130) 


and finally from equation (114) comes D t (r f /R Q ). The slopes given by equations (128) and 
(123) are matched (eq. (127)); and the location of the freezing interface r f /R Q , the time 
Tp and the surface temperature T(0) are now determined. An example appears in a 
appendix F. 
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Of particular interest is the time required for the surface temperature to reach the 
Leidenfrost temperature T^. From equation (126), 


t L 



T L ” T s 


T(0) - T c 


(131) 


Equation (131), while limited by simplifying assumptions, gives an estimate of how long 
the sphere will float. 


1 
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APPENDIX J 


SAMPLE CALCULATIONS 

In order to test the theoretical model, an experiment was performed using water and 
liquid nitrogen. The problem considered here is to calculate the observed cnanges in the 
characteristics of the sphere. 

Consider a sphere, radius R q , of dye-tinted water "floating” on a sea of liquid ni- 
trogen. As the fluid within the sphere freezes, the color changes; it is estimated that at 
least a 0. 1R Q layer of ice is required to make the color change vivid. Because the 
sphere of water must contact the nitrogen to generate its vapor cushion, the surface 
should become crystalline; that is, ice will form. The crystalline structure forms a non- 
uniform surface; ice spires protrude into the encapsulated liquid and into the vapor gap, 
thereby initiating an early destruction of the vapor cushion. This description serves only 
to remind the reader that these calculations are greatly simplified with respect to the 
physics of the problem. With these limitations in mind, consider an 0. 305-centimeter 
sphere of water floating on liquid nitrogen, which possesses a 6* of approximately 
113-j^ , as is shown in appendix H and illustrated in figure 10 for R q = 0. 3 centimeter. 
Therefore, the heat-transfer coefficient will be calculated from equation (C25), 9* < v. 

The water temperature at inception is 298 K, and the saturation temperature of 
liquid nitrogen is 77.4 K. The drop surface temperature is assumed to be 273 K. There- 
fore, the mean drop temperature is 286 K, and the vapor properties are evaluated at 
175 K. These properties, obtained from references 7 to 9, are given in table IV. 

Calculating the parameters W d , w*, and S in equation (C25) yields 

W d = — p d R^ JL = i2L (980)(0. 0284) = 116. 5 dynes (Jl) 

1 3 g c 2 3 * 

2 2 

W* = r2 = (0. 00196) f— 16x10 ^2 = 6. 86xlO" 6 dynes (J2) 

g c ° \ 0.00196 ) r 2 


+ I( C P iT ) 

= 198-u +2 

Rl . 045)(196)i| 

A x J 

l 2 

L 198 Jj 



= 307 Joules/cm 


1.031 


(J3) 

(J4) 
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TABLE IV. - PROPERTIES OF VAPOR, WATER SPHERE, 


AND LIQUID NITROGEN 


Vapor properties: 


Viscosity, /i, dyne-sec/cm 

1, 16x10"^ 

Thermal conductivity, k, W/(cm)(K) 

1. 76xl0" 4 

Density, p, g/cc 

0.00196 

Specific heat, C^, J/(g)(K) 

1.045 

Properties of sphere: 


Thermal conductivity, W/(cm)(K): 


Liquid, 
Ice ’ k ice 

0.0057 

0.0235 

Density, g/cc: 


Sphere, p d 

1 

Ice ’ Pice 

Specific heat, J/(g)(K): 

0.9 

Sphere, (c p ) 

4. 19 

Ice ’ ( C p)- 
' w ice 

1.9 

Latent heat of fusion, y, J/g 

336 

Radius of sphere, R q , cm 
Properties of liquid nitrogen: 

0.305 

Density, p^ , g/cc 

Latent heat of vaporization, J/g 

0.81 

198 

Acceleration of local gravity, g- cm/sec 2 

980 

Gravitational constant, g c , g-cm/dyne-sec^ 

1 


9 lO I® 

f( 0 *) = cos ' 2 6* - 2 cos e* + 1 = cos 66 -^ + 2 cos 662 + 1 = 1. 956 


(J5) 


Using these values, the Nusselt number may be calculated according to equation (C25) as 


h ,R 

N Uj = -1^ = 1 + 

d 


= 1 + 


W, 


677W*Sf(0*) 


1/4 


116. 5 


677 ( 6 . 86x10 -6 )(1. 031)(1. 956) 


1/4 


= 26. 8 


(J6) 


The heat-transfer coefficient is readily determined 
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h d = — Nu d = ( 1 - - ? 6Xl °-\ 26. 8) = 0.0155 watt/(cm 2 )(K) (J7) 

R Q \ 0.305 / 

The average heat-transfer coefficient to a 0.6-centimeter-diameter sphere warming 
in a gas stream (ambient air) flowing at velocities to 1/4 meter per second was found 
experimentally to be approximately 


Nu amb = 10 


(J8) 


Therefore, for these calculations it was assumed that 


h amb = — Nu amb =- 76X10 4 (10) = 0. 0058 watt/(cm 2 )(K) 
R Q mD 0.305 


(J9) 


The average heat transfer for the floating sphere would be dependent on the surface area 
exposed to gaseous nitrogen A amb and that exposed to the Leidenfrost phenomenon 


h = 



, ^amb , 
+ h 


amb 


(J10) 


For the case of the 0. 3-centimeter-diameter sphere, discussed in appendix H, where 
9 * = 113-^ , the heat-transfer areas become 

A fb = -2 ttR 2 (cos 9 * - 1) = 4ttR 2 (0. 7) (Jll) 

A amb = 1 - A fb - 4 * R o<°- 3 > ( J12 > 

These areas are not significantly altered for the 0. 305- centimeter-diameter sphere. 

Substituting equations (Jll), (J12), (J9), and (J7) into equation (J10) yields the aver- 
age heat-transfer coefficient for the floating sphere 

h = (0. 7)(0. 0155) + (0.3)(0 0058) = 0.0126 watt/(cm 2 )(K) (J13) 


The time required to freeze a spherical shell 0. 1 R q thick may now be determined 
by using equation (121). Evaluating the constants 1 C and 2 C gives 
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i c _ i _ R o h d = j _ (0. 305)(0. 0126) = Q g37 


k. 


ice 


0.0235 


(J14) 


o R^h , 

Z C = - 0 Q = 0. 163 

k. 

ice 


(J15) 


Before u can be evaluated, the heat-transfer coefficient h^ must be selected. In this 
part of the calculation, assume it to be k^/R Q . 


- 3 C = h ^ d ' Tf) = ** (Td ' Tf> = 0.0057(286 - 273) = ^ 002634 se(; -l 


p ice yR o p ice r R o 0.9(336)(0.305) 2 


(J16) 


^ _ h d (T f - T s } _ 0.0126(273 - 77) _ Q 027 sec -l 
P ice yR 0 0. 9(336)(0. 305) 


(J17) 


2 

C R = ( 2<:3c ) - 4 WC = 0. 163(-0. 00263) 2 - 4(0. 84)(-0. 00263)(0. 027) 


= 2. 36xl0' 4 sec" 2 


(J 18 ) 


As the characteristic C R > 0, equation (12 la) is applicable, or 


t = -- 


1 


0.00263 


(1 _ 0 9) i °- 027 ln | (l-8)(0 837)(-0. 00263) + (0. 163)(-0. 00263) - 0.0154 
. " ‘ + 0 . 0 1 54 1(1. 8)(0. 837)(-0. 00263) + (0. 163)(-0. 00263) + 0.0154 


x 2(0. 837)(-0, 00263) + (0. 163)(-0. 00263) +0. 0154 1" 
2(0. 837)(-0. 00263) + (0. 163)(-0. 00263) - 0. 0154|_ 


= 3.7 sec 


(J19) 


This is the time required for the entire sphere with a liquid core heat source to freeze 


a 0. 1 R q shell. 


Assume now that the sensible heat of the liquid core is negligible; the time as calcu- 
lated by equation (122) becomes 
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t - 0- 9(336)(0. 305) U _ (Q g) 3 ^ 0.305(0,0126) [j 3(0 g) 2 
3(0. 0126)(196) l ' 2(0.0235) L ' 



= 3.2 sec (J20) 

The difference in time with and without sensible heat is about 0. 5 second, about a 
15 percent difference. The time required to freeze other thicknesses by these techniques 
is shown in figure 12. 



Figure 12. - Freezing interface advancement time as function of position. 


The time required to obtain an 0. 1 R q shell of ice may also be obtained through the 
simplified energy balance. If we let 


^bs ^sen + ^ce 


(J21) 


then the time to remove the sensible heat without any surface freezing is given by 
equation (32). 
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( 32 ) 


mCj (T h 0 - T f ) 

t V P/d d >° f _ (0.117)(4. 19)(298 - 273) _ 4 2 see 

Sen (hA) d (T f - T g ) [(0.0126)(4 jt)(0. 305)^(273 - 77) 

The time required to freeze the 0. 1R Q thickness shell is given by equation (33): 


R_p { 

Lml 

r 

cr ice 

L \ R J _ 

(0. 305)(0. 9) [l - (0.9)^1(336) 


ice (hA) d (T f - T g ) 3h d (T f - T s ) 3(0. 0126)(273 - 77) 


= 3.4 sec 


(33) 


Based on these estimates, the time required to freeze a 0. 1R Q concentric shell would 
be (eq. (J21)) 


obs 


= t 


sen 


+ t ice = 7 - 6sec 


(J22) 


Heat-transfer 
coefficient to sphere, 



Figure 13. - Solution for changeover point from equation (121) to equation (126). 
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In comparing this time with equation (J19), it is evident that all the sensible heat cannot 
be extracted without freezing at the surface. These techniques are compared in figure 12 
at various positions of the freezing interface. 

In figure 12, it may be noted that the calculations were carried to r^/R^ = 0 without 
regard to the changeover point as discussed in appendix I. Figure 13 illustrates the vari- 
ation in the temperature-time gradients D^.(T _ T ) for equations (121) and (126) as func- 
tions of time and surface temperature. The changeover occurs where the slopes are 
equal: 


D. (t - T ) = D t (T - T a ) 

tV s 'eq.(I21) ^ s 'eq.(I26) 

The required parameters T(0) and are then determined from figure 13. 

T(0) = 238 K 
Tj = 12. 9 sec 


(127) 


(J23) 


To estimate the time required for the surface temperature to drop to the ordinary 
Leidenfrost point (i. e. , the termination of the floating conditions relative to a smooth 
steel surface), equations (131) and (J23) are used, with T^ = 98 K, 


t L = " 


3h,, 


In 


T L" T s 


T(0) - T, 


(0. 305)(0. 9)(1. 9) ln / 20 \ 
(3)(0. 0126) \160. 6/ 


= 28.3 sec (J24) 

The time from placement of the liquid-water sphere on the liquid nitrogen to the time 
it will sink is, from equations (124), (J23), and (J24), 

t L = 12.9 + 28.3 = 41.2 sec (J25) 

The result is in fair agreement with the data of table III. 

The pseudo -steady- state analysis and the simplified energy balance techniques vary 
in their agreement with the data. The effects of things we do not know, such as the 
Leidenfrost temperature for a water sphere on liquid nitrogen, and the simplifying as- 
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sumptions of appendix I indicate (1) the pseudo-steady-state analysis should be modified, 
and (2) that the experimental data should include thermocoupled spheres to determine the 
Leidenfrost temperature and the degree of subcooling of the sphere during freezing. 

The calculations made in this appendix are summarized in table V. 


TABLE V. - SUMMARY OF SAMPLE CALCULATIONS 


Parameter 

Calculated value 

Experimental value 

Sphere weight, W d , dynes 

116. 5 


Weight for nondimensionalizing, 

6. 87xl0' 6 

- 

w*, dynes 



Sensible-to- latent-heat parameter, S 

1.031 


Nusselt number, Nu 

26.8 


Heat-transfer coefficient to sphere, 

0.0126 


h d , W/(cm 2 )(K) 



Pseudo-steady-state solution 


Time required to freeze a 0. 1 R q 
concentric shell, t, sec: 



With sensible heat 

3.7 


Without sensible heat 

3.2 

8 to 9 

Time sphere will float, t^, sec 

41.2 

23 to 27 

Simplified energy balance 


Time require to freeze a base 0. 1 R q 

7.6 

8 to 9 

thick, extracting all sensible heat 
first ’ ‘obs- sec 



Time sphere will float, t^, sec 

46.0 

23 to 27 
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APPENDIX K 


ESTIMATE OF LEIDENFROST TEMPERATURE 

The Leidenfrost temperature is the temperature on the pool-boiling curve associated 
with the change from film to transition boiling. A recent, but as yet unpublished, work 
by Baumeister and Henry shows how large changes in the surface characteristics affect 
the measurement of the Leidenfrost temperature. They show that the Leidenfrost tem- 
perature of saturated water on glass is considerably higher than the Leidenfrost tempera- 
ture of water on steel. These results are shown in figure 14 in terms of the temperature 
differences above the liquid saturation temperature . 

For liquid nitrogen on a metallic surface, the Leidenfrost temperature is about 97 K, 
or about 20 K over the saturation temperature of nitrogen at atmospheric pressure. In 
this analysis, however, the Leidenfrost temperature of an ice-nitrogen combination is 
required. At present, there are no available data from which to estimate this Leidenfrost 
temperature. Consequently, an estimate of the ice-nitrogen Leidenfrost point is made 
from the data shown in figure 14. 

First, the dimensionless Leidenfrost temperature difference 0^ for smooth or 
rough surfaces, shown in figure 14, is assumed to vary linearly as a function of the nat- 
ural logarithm £/£ s teel‘ However, the variation of £ with surface roughness, a geom- 
etry factor, has not yet been assessed. Consequently, the superposition principle will 
be used herein. For smooth-surfaced ice, the ratio of ^ice^steel * s *3.7. Conse- 
quently, the Leidenfrost temperature of smooth-surfaced ice on nitrogen will be approxi- 
mately 1. 89 times the value measured on steel. 



Figure 14. - Effect of surface on Leidenfrost temperature. 
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Besides the correction for surface properties, an additional correction for surface 
roughness must be considered in estimating the Leidenfrost point of an ice-nitrogen com- 
bination. The previous estimate for the Leidenfrost point assumed a smooth surface on 
the ice sphere. However, as was mentioned in the body of this report, ice crystal pro- 
trusions occur on the surface. Cumo, et al. (ref. 13) investigated the effect of surface 
roughness on the Leidenfrost point. Cumo, et al. (ref. 13) report that the Leidenfrost 
temperature difference for water on a smooth gold-plated surface is 285° C, while the 
Leidenfrost temperature difference on a sandblasted surface is 335° C. Thus, the ratio 
of the temperature differences for a rough to a smooth surface is 1. 27. 

The increase in the Leidenfrost temperature difference is a result of the decrease 
in the effective value of £ at the surface due to the surface protrusions that penetrate 
through the vapor layer. At present, however, there is no means of estimating the new 
value of Consequently, it is assumed herein that the dimensionless Leidenfrost dif- 
ference for a rough surface will be proportional to that measured on a smooth surface. 
The constant of proportionality is taken from the data of reference 13 to be 1.27. This 
constant is assumed to be independent of the value of |. Multiplying the dimensionless 
Leidenfrost temperature difference for the smooth surface by this value gives the upper 
line in figure 14 which applies for rough surfaces. 

From figure 14, Leidenfrost temperature difference for the rough ice-nitrogen com- 
bination is now estimated to be 2.4 times the value measured on steel, about 126 K. 

It must be emphasized that this is only a first-order estimate. The data and line 
shown in figure 14 do not take into account property variations of the liquid, only the 
metallic surface properties are considered. 


« 
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Motion-picture film supplement C-267 is available on loan. Requests will be filled 
in the order received. 

The film (16mm, 14 min, color, sound) presents background material and discusses 
the forces acting to levitate a sphere in Leidenfrost boiling on a liquid interface. The 
associated phenomena are demonstrated by using blue-tinted water spheres floating on a 
sea of liquid nitrogen. As the sphere freezes, the color changes from blue to lavender 
and the partly frozen sphere is free to roll about on the surface because of the shift in 
the center of gravity. The freezing sphere crystallizes and frequently splits into two 
major fragments. These fragments are at first repelled and then attracted until they 
reattach. As time progresses, the surface temperature decreases until transition boil- 
ing begins; this destroys the vapor cushion and the sphere sinks to the bottom. Water 
spheres can be dropped several centimeters and still float; however, size, fluid prop- 
erties, and dropping distance are important. A sequence of two floating dye-tinted spheres 
which attach is followed to the point where they sink. 

The phenomena are again demonstrated using n-heptane. A small-scale fuel spill is 
simulated. 

The whole process is then reversed (i.e. , the liquid nitrogen is floated on water). 
Vapor flow patterns are very easily detected and ice paddy residues are common. The 
liquid-nitrogen spheres are seen to vibrate in distinct modes, and ice tracers enable the 
observer to follow flows internal to the liquid- nitrogen sphere. 

The process is discussed as a flow visualization technique. 

Film supplement C-267 is available on request to 

Chief, Management Services Division (5-5) 

National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, Ohio 44135 
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Washington, D.C. 20546 
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